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ABSTRACT 


Title of Thesis: Relativistic Electron Precipitation in the Auroral 

Zone 

David Simons, Doctor of Philosophy, 1975 
Thesis directed by: Professor David L. Matthews 

The energy spectra and pitch angle distributions of electrons in the 
energy range 50 keV to 2 MeV have been determined by a solid state 
electron energy spectrometer during the Relativistic Electron Precipi- 
tation (REP) event of 31 May 1972. The experiment was carried aboard a 
Mike-Cajun sounding rocket as the University of Maryland component of a 
joint American-Norwegian (WASA-HDRE) ionospheric investigation. 

The difficulty of determining the expected electron flux prior to 
the experiment required an instrument with a large dynamic range. The 
design and theoretical modeling of this instrument is described in 
great detail. The extensive preflight calibrations on the NASA-Goddard 
2 MeV Van de Graaff electron-proton accelerator are compared with the 
theoretical model. 

The electron pitch single distributions are determined from a know- 
ledge of the rocket aspect and the direction in space of the Earth's 
magnetic field. The rocket aspect determination is therefore treated in 
depth and a method is developed to compensate for the malfunctioning 
of the aspect magnetometer. 

The electron fluxes during the REP event were highly variable demon- 
strating correlated energy, flux and pitch angle pulsations with time 
periods less than one second. Increases in flux wei*e accompanied by marked 
filling of the loss cone at lovrer energies (near 50 keV). Within the time 




resolution of the experiment, these fliixes were simultaneous with a 
w hardening of the energy spectrum characterized by marked increases of 

the near relativistic electrons. The higher energy particles showed a 
much more pronounced loss cone than the lower energy particles during 
■ these flux increases. The procedures for determining these results 

i from raw data count to final conclusion are also described. 

I 

i Drawing upon the q.uasilinear equations of plasma wave-electron 

i interactions, a theoretical model for the production of relativistic 

I electrons is proposed. A self consistent set of fully relativistic 

: equations for the evolution of the electron distribution function due 

i 

to the interaction of the electrons with parallel propagating whistler 

j waves is derived in the Appendix. An examination of these equations leads 

I 

} to the conclusion that at comparatively low background electron densities, 

I 

I . 

f the anomalous doppler resonance leads to the acceleration of near rela- 

I 

I tivistic particles. The results of a computer solution of the five 

i coupled integrodifferential quasilinear equations confirms this conclusion. 

I 
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CHAPTER I 


INTRODUCTION 


The REP Event 

The Relativistic Electron Precipitation (rEP) event is a subclass 
of daytime auroral absorption belonging to the Slowly Varying Absorption 
(SVA) class of absorption events. It has a number of features that have 
set it aside as a fairly special form of magnetospheric-ionospheric 
phenomenon. Auroral absorption is the abnormal absorption of cosmic 
ratio noise (CRN) in the ionosphere of the auroral zone as determined by 
a riometer (Relative Ionospheric Opacity meter) . Absorption is caused 
by increased ionization in the lower ionosphere due to x-rays and pre- 
cipitating energetic charged pairtioles. The increased ionization leads 
to resistive dis'-.ipation of the transmitted signal (i.e., ionospheric 
absorption). An absorption event is any period of marked increased 
absozption . 

Ever since Little and Leinbach (1959) developed a sufficiently 
stable receiver to make absoiption measurements over a long period mean- 
ingful, ionospheric observers have been using the riometer to monitor 
ionospheric activity. , Absorption events are normally characterized by 
the physical appearance of the printed riometer record. Hence the name 
slowly varying absorption characterizes the slow smooth decrease in 
signal intensity extending over as much as an hour and the tisually slower 
recoveiy (one or two hours) to the normal level. This slowness during a 
daytime event is due in part to the sluggishness of the daytime iono- 
sphere. As pointed out by Ansari (Ansari* 19^5) this behavior sets these 
daytime events quite apart from the typical ni^ttime auroral absorption 


vhich is characterized hy rapid fluctuations in the signal intensity. 

The daytime SVA has heen discussed extensively oy Ansari (1965)9 Ortner 
and Riedler (1964) and Bewersdorff et al. (1968). 

The REP event was first distinguished from other SVA's hy Bailey 
and Pomsrantz (1965) using a riometer in association with VHF forward 
scatter radar. The forward scatter experiment gives some idea of the 
depth of penetration into the atmosphere of the ionizing radiation 
which is causing the absorption. 

Bailey and Pomerantz determined that ionizing particles were pene- 
trating into the atmosphere to a depth below TO km during certain SVA's. 
The fluxes of charged particles below the reflection hei^t of the VHF 
forward scatter radar were sxifficient to cause very large absolution of 
the scattered signal. By modeling the absorption mechanism they inferred 
that the effective atmospheric cutoff energy for electrons causing this 
absorption was about UOO Kev. Protons causing the same absorption would 
have had a cutoff of 10 Mev and were therefore considered much less 
likely. Bailey and Pomerantz therefore coined the expression Relativ- 
istic Electron Precipitation. 

The REP is a distinctly daytime event meaning that the geograph- 
ically widespread absorption is observed in the dayside auroral iono- 
sphere. There is generally little or no local magnetic activity 
observed, (Ansari, 1965s Bailey, 1968) and REP's occur following the 
ea^ansive phase of magnetosplieric substorms as observed on the ni^t- 
side of the Earth (Thorne, 1974). In order to put this in perspective 
we must dwell on the structure cf the magnetosphere and define a few 


terms 
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I.B. 'A Brief Reviev of MaRnetospheric Structure« 

The solar plasma streaming by the Earth encloses the Earth's 
magnetic field in a cavity vhich ■^’^e call the Magnetosphere. This essen- 
tially fluid type interaction of plasma and field draws out the Earth’s 

* 

field on the ni^tside into a long tail ari compresses the field on the 

dayside confining it to something of the order of 10-12 Earth radii (Rgi 

Within about 6 R„ the magnetic field bears some resemblance to a dipole 

confi, juration (Dessler, I968) . Farther out than 6 Rg field lines are 

extended toward the tail on the ni^tside and compressed on the dayside. 

A magnetic neutral sheet extends up the tail from about 10 Eg and has 

been observed to exist well beyond 80 (Ness* 1969). Field lines that 

extend into the tail past the beginning of the neutral sheet are called 

open field lines and may very well merge with the solar wind magnetic 

field some point far downstream. It has not been established whether 

particles from the solar ward have direct access into the tail. The 

neutral sheet and the central portion of the tail are hathed in a plasma 

referred to as the plasma sheet. Within aborrt 4 Rg of the Earth is the 

stable trapping region which contains a plasma which is shielded from 

the fairly extreme fluctuations which are observed in association with 

magnet ospheric storms, {Van Allen, I969), This region has a cold plasma 

• • 3 

density of about 100/ cm . Except during storm conditions, there is a 
precipitous drop off in density between Ej, and 5 Eg called the 

q 

plasmapause (Carpenter, I966) where the density falls off to about 1/cm . 
The region between 5 Eg ^d the plasma sheet is the region of the 
magnetosphere that is connected hy the magnetic field to the ionosphere 
of the A^lroral Zone, Tliis is a transition region between closed field 
lines and field lines that extend on the nightside into the magnetotail. 


k 

Drastic changes in the magnetic field and particle population of this re- 
gion acconrpany magnetic snhstornis (see helow) keeping the region in a 
state of change. Even during relatively quiet periods there seems to be 
continual acceleration of particles in this region £is manifested by 
luminosity in the auroral zone. The auroral activity is caused by par- 
ticles precipitating into the ionosphere at the feet of the field lines 
(Akasofu, 1968 ). 

Particles trapped in the Earth's magnetic field drift aroxmd the 
Earth due to the gradient and curvature in the magnetic field (gradient 
drift). These drifting particles effectively connect opposite sides of 
the magnetosphere into regions called drift shells. Barring perturbations 
that distrirh the adiabatic motion of the particles, a particle on the night- 
side of the matnetosphere should drift around the Earth and return to 
the nightisde remaining on the same drift shell. In this manner magnetic 
field lines that penetrate the transition region and even extend out into 
the magnetosphere have effective communication (by drift shell particles) 
with the field lines on the dayside. The unsteady state nature of the 
region external to the nightside plasmapause is carried around to the 
dayside. The field lines from this region define an area in the Earth's 
ionosphere called the Auroral Oval (Feldstein, 1969^ Akasofu, 1968). 

During times of solar activity the solar wind carries disturbances 
to the magnetosphere. Particularly large magnetospheric disturbances 
are called magnetospheric storms. Flare- associated energetic particles 
reaching the Earth's magnetosphere are associated with marked compres- 
sions of the magnetic field. These compressions are often followed by 
successive periods of marked magnetic activity which last, for times of 
one to three hours. These one to three hour events are called substorms. 
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There are usually several suhstorros associated with each magnetospheric 
storm (Ahasofu, I 968 ) . As determined hy ntanerous satellite observations 
there are marked increases in the trapped energetic particle populations 
(E > 1 Kev) of electrons and protons and a corresponding increase in the 
am’oral activity in the Oval. Large electric fields are induced in the 
auroral ionosphere causing concentrated currents in the region of the 
Oval. These cxnrrents are called the atiroral electrojet. The currents 
are detected on the groiuid by the severe fluctuations they cause in the 
local magnetic field. Ground based magnetometers which monitor the 
vector components of the local magnetic field record changes of several 
hundred Y (lY = 10*”^ gauss) in the different components of the field. 
The most intense precipitation and the associated magnetic activity are 
generally observed in the midnight sector of the Aiiroral zone. As 
determined by rocket in-situ observations of electrons and protons. 

X-ray observations made from high atmosphere balloons, and satellite 
particle observations, the majority of cosmic radio noise abso3T>tion in 
the ionosphere is contributed by electrons in the 10 Kev energy range. 
There are significant increases in the trapped electron population for 
particles in excess of 50 Kev (Lezniak et al. I 968 , Lezniak and Winckler, 
1970 ) as well as protons in the 1 to 50 Kev energy range although these 
particles do not make a major contribution to auroral absorption. It 
has been suggested by Lezniak and Winckler (1970) that these energetic 
particles are injected from the tail do^m into the trapping region near 
local midnight and then drift from there to other local tiroes. 


• 4 . 



With this "brief introduction we can atten 5 >t to place the REP in some I;; 

context as a s\ 'jstorm-related phenomenon occurring in the morning side 
ionosphere. The corresponding magnetospheric region would presunably f? 

• r 

receive excess trapped particles from gradient -B drifting particle 



populations injected on the night side near midni^t and not directly 
from the plasma sheet. Little or no magnetic activity suggests there 
are no local current systems and little, if any, distortion in the local 
magnetic field. The miisiial hardness of the energy spectrum may indi- 
cate that the precipitation mechanism for REP*s is significantly different 
from that responsible for nights ide precipitation. Looking at Pigure 1 
of the auroral zone and auroral oval, we see that dayside auroral-zone 
latitudes actually can "be below the auroral oval and are connected there- 
fore to field lines that are closed in the magnetosphere. The REP is 
reported to extend over a region as great as 1,000 Kin (Thorne, 1974) in 
extent and therefore involves a significant energy source. 

The highly unusual nature of a REP in terms of other auroral-zone 
events makes it a very interesting phenomenon to the auroral and magneto- 
spheric physicist. The possibility that a detailed study of the phe- 
nomenon might shed some new light on the entire field of complex 

* 

ionospheric-magnetospheric interactions is a veiy appealing one. Ground 
based measurements such as the VHP forward scatter and riometer only 
give qualitative indications as to the nature of the precipitating flux 
while balloon measurements of the Bremstrohlung X-rays can say nothing 
about the directional distribution of the particle fluxes about the 
magnetic field, A properly designed in-situ rocket measurement of the 
actual particle fluxes resolves the question of particle species, energy 



U. 


Figure 1. The Auroral Oval 

A projection of the auroral oval onto the Earth's 
surface at about 6:00 UT taken from Akasofu (1968), 
Northern Scandinavia is in the nighttime auroral oval 
but may be seen to be well south of it according to this 
presentation in the morning hours. 
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dependence, and pitch angle distribntions, and allows for high reso- 
lution time analysis. It was hoped that such a rochet-hased REP escperi 
ment could yield some insight into the complex plasma physical processes 
which energize thr particles that precipitate at various times and places 
into the uppex* atmosphere. 
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I.D, The Orpianization of This Reports 

o 

The following detailed discvission of the EEP experiment is divided 
into four chapters. Chapters 1 and 2 deal with the experimental con- 
siderations. Chapter 1 describes the design, construction and calibra- 
tion of the electron spectrometers while Chapter 2 details the solution 
of the rochet attitude problem. For those interested principally in 
the results of the experiments these chapters may be shipped. Chapter 3 
describes the data obtained during the EEP of May 31, 1972, and the 
reduction of that data. In Chapter U we disctiss these experimental re- 
sults in the li^t of present magnetospheric theories. 
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CHAPTER II 

THE DESIGN AND CALIBRATION OF A NINE 
CHANNEL ELECTRON ENERGY SPECTROMETER 

II.A, Choosing An Electron Detector 

1 . Crii^eria For Counting Capahility 

Bail^ (1968) wrote an extensive review of the REP based upon 
three years (1964-I966) of ground "based observations in Alaska using the 
riometer and forward scatter experiments. Bailey's results provided a 
major source of particle flux data used in establishing design criteria 
for our in-situ REP experiment. In view of this we will discuss this 
data before discussing the design of the spectrometer. 

The differentiaJ. absorption of cosmic radio noise at a given 
hei^t in the ionosphere is a function of the local ionization and the 
local collision rate between electrons and neutral particles (Sen and 
NyHer, 1960) . Even if one assumes a fixed collision rate height pro- 
file, various ionizati'on hei^t profiles may produce the same total 
absorption (at a sin^e frequency) observed on the ground. By adding 
the foiward scatter data to the riometer data the ambiguity may be re- 
moved provided the scattering height is well known. Bailey could'not 
fix the scattering height, but by assuming vario^ls heights he could 
produce a range of flux and hardness parameters for the REP's that he 
observed. The flux of precipitating particles was assumed to be iso- 
tropic over the downward hemisphere with a unidirectional integral flux 
of the form 

-E/E 

J(>E) = J e ° 
o 



I 


\ 

I 

H 

H 


2*1 
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(Jj with J heing the number of particles per cm^ per se:. per star, above 

E (in Kev) and E (called the e-folding energy in 1st) being the 

o 

spectral hardness parameter. In a very crude sense, T, controls the 
amount of ionizaticaj produced in the ionosphere and Z. controls the 
depth of ionization production. 

As an exan^le of Bailey’s results we consider the REP of 
November 13, 1965* Choosing a scattering height of JC fan for the for- 
ward scatter produced a spectra with E^ = 68 Kev and - 7*6 x 
5 2 

10'^/(cm X sec. x ster. ), while choosing a scattering height of 75 Km 
produced a spectrum with = 34 Kev and “ 4.8 1 10^(cm^ sec ster). 
So appreciate the difficulty this presented us with Then designing the 
spectrometer, let ijs look at the flux at several energies. For the 
softer of the two spectra the flux at 30 Kev is four times the flux at 
30 Kev for the harder spectrum. At 400 Kev however, the flxcc of the 
harder spectrum is some 56 times the flux of the softer spectrum. All of 
the REP’s analyzed by the forward scatter-riometer technique 
same wide range of possible flux values. We were therefore forced to 
design with the widest possible dynamic range in mind. The major con- 
sideration for oin: device was its abilii^ to cope with as large a range 
of hardness and flux parameters as possible. 

2. The Advantages Of Solid State Detectors 

Charged particle detection on a sounding rocket calls for com- 
paratively small detection devices with high reliability. The semi- 
/ fj conductor type particle detector was the natural choice in the 50 Kev 

to several Mev range. Such detectors can he small and have straight 
forward geometry, while energy resolution in this range tends to be 
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limited more "by space available for associated electronics and telemetry 
i-J' information rate, than by the device itself. We wished to be able to 

detect electrons with energies less than about 2 Mev and greater than 
some 50 Kev and in addition had to be prepared for a wide dynamic range 
of particle fluxes . 

The semi-conductor detector is essentially a slab of silicon 

' a 

(or germanium) in which a charged particle stops through ionizing inter- 
actions resulting in the production of mobile charge carriers in the 
conduction band of the silicon. These charges are then collected by an 
electric field which is maintained across the detector. Superficially 
this is veiy much the same behavior as a gas propojrtional counter. The 
theory of operation is fairly well developed and the literature is ex- 
tensive. (see Taylor j I963 or Deamaley and Northrop, I966), We will 
treat the basic operation of the detector only veiy briefly. 

Several different types of semi-conductor detectors are made. 
Our choice of a lithium compensated (P.I.N.) detector was dictated by 
the con^aratively large depth of silicon necessary to stop 2 Mev 
electrons and because of the need of an energy threshold below 50 Kev. 

The average integrated path length of an electron coming to 
rest in a material may be determined by use of the Bethe stopping power 
formula. (Bethe, I938), 


dX 


1 U 
Uttb ZN 


B 


mv 


( 2 . 2 ) 



with B given by 



2 

mv E 

2 I^( 1 -B^) 


- I In 2 - 1 + B^) + 1 - B^ 


( 2 . 3 ) 
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This formula is good for energies below 10 Mev where ionization is the 
principal energy loss mechanisin. Z is the atomic mimber of the absorber, 

3 

N the mmber of absorber atoms per cm , m the mass of the electron 
and I the mean ionization potential of the absorber. This formula 
yields a path length which is about 1^0^ greater than the average depth 
of penetration because of the considerable Isirge-angle scattering. Us- 
ing the appropriate values for silicon it is easily determined that 
somewhat over two mm are required to stop one-Mev electrons, 

3. Some Properties f Silicon Detectors 

In barrier layer detectors of the diffused junction and sur- 
face barrier type the depletion region is obtained by applying a large 
voltage. The larger the voltage the deeper the depletion region will be. 
The hipest resistivity silicon (i.e. the purest and most expensive) will 
tolerate fields of a thousand volts/mm at best before breaMng do^m 
completely. Depletion depths of a few millimeters have been achieved 
but at the expense of large leakage curi'ents because of the large re- 
verse bias required. A far more inexpensive device which operates at 
much lovrer voltages is the lithium compensated silicon detector. The 
lower voltages mean less leakage current and therefore less noise, 

P.I.H. is an abreviation for p-type-intrinsic-N-type semi- 
conductor detector and describes the physical arrangement of the three 
basic regions in the lithium compensated silicon detector. A slab of 
XHtype silicon crystal is painted on one surface with an oil suspension 
of lithium. The lithium is then diffused under high teirperature into 
the crystal. lithium, a donor inpurity, greatly overcompensates the 

p-type acceptors at one surface of the crystal and results in the forai- 


ation of an N~P Junction at about 300 microns in the crystal. 

Temporaiy electrodes are plated to the ciystal surfaces and a reverse 
bias is then applied across the Junction 30C volts). The lithium 
drifts all the way across the silicon creating a large region of silicon 
with almost the properties of pure and undoped silicon. The crystal now 
has an n-type region some 200 to 300 microns deep and an effective 
intrinsic region with no majority carriers the rest of the way across. 

The surfaces are etched ^wn and a surface Junction is obtained on the 
intrinsic surface by the application of a thin layer of gold. 

(30 ugm/cm ). Kiis surface is the entrance surface for charged particles 
and constitutes a small dead layer. Our detectors had about l60 ^ 

(31 ugn/cm ) dead layer of gold. The other surface after etching is 
coated with aluminum. The intrinsic region (in our detectors about 3 
mm thick) is che active part of the detector. In Figure 2.1 we show 
the detector in its motmting. 
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H.B. The Collimator 

The solid state detector vas enclosed in a Brass and lead coll- 
mator. The collimator geometry was determined "by considering the 
possible intensities of electron precipitation that mi^t occur during 
rep's as veil as the angular resolution set hy the telemetry and rocket 
spin. 

1. Geometrical Factor 

The differential flux of particles in a given direction, 

F{u,x,E) , sometimes called the unidirectional f lux or directional 

_2 -2 -1 -1 

intensity, is usually expressed in units of cm star Kev sec 
X is a position vector in space, u is the unit vector in the direc- 
tion of the particle motion, E is the particle energy, and F(u,x,E) 
points along u . The counting rate N , in units of sec”^ , of a 
collimator-detector combination is a function of the geometry of the unit 
and the efficiency of the detector. For the collimator-detector 
geometry shown in Figure 2.2, E is given by 


N = 


"^2 

mi 


(?*da^){Ti*da2) E(e) 


p 

DtD” 


dEdt 


( 2 . 4 ) 







tiiilirY-imil 




. 18 


vhere the geometrical terms are defined in Figure 2.2 E(E) is the 
detector efficiency, which ve will treat fully in section IC.D. life is 
the counting interval. Since the collimator does not distinguish he- 
tween directions in the field of view we ntust assume some average fltix 
over the aperature. Letting 


F(E) = <'|f(u,x,E) j)>^ ^ 

x,n,t 

we may write equation (2.4) 


H = G 


I F(B) E(E) 

E 


dE 


where G , called the geometrical factor, is strictly a function of the 
collimator-detector geometry. Using polar coordinates r^ , and 
*2 * *^2 planes of and Ag (see Figure 2.2) separated by 

distance h we find that G may be expressed as 


2w 27T 2„ 

fl f2 cos ar,r. 


f r fl cos ar^r 

I 1 i 


0 0 00 0 


o 2 2 2 li 

with D = r^ + Tg - 2r^rg cos((j)^-(^g) + h and cos a = — . Stflj- 

stituting in these values and performing the integration yields 




2 2 
+B|+h‘" 


- r, 


/_2.„2.. 2.2 

(Rj+Eg+h ) 


**E^Rg 


} 


(2.5) 


Since it is very difficult to produce an isotropic flux of 
energetic electrons in the laboratory, only indirect means of cheeking 
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the geometriceJ. factor nay he resonably lised. This point is discussed 
fully in section II ,D, 

The value that one chooses for the geometrical factor is deter- 
mined hy the expected flux levels, accumulation time interval, the width 
of the energy channels of the energy analyzer and the capacity of the 
data accumulators. Following Bailey (1968) we modeled the expected 
differential energy spectra in the form 

J 

F = “ exp {-E/E ) (2.6) 

Q 

vith vaiying from 25 Ker to 150 Kev and varying from 10 to 
6 *2 —1 *"1 

10 (cmT sec ster” ) - Obviously one cannot adequately cover all these 
possibilities with a single device and achieve optimum performance for 
all possibilities. We designed for optimum conditions toward the center 
of the ranges and attenpted to extend our limit as far as possible in 
both directions on hardness and flux levels. The hardest spectrum here 
yields about 15/ (cm x sec x ster x Kev) at 1.0 Mev while the softest 
spectrum gives essentially nothing at 1,0 Mev. We planned for an inte- 
gral channel as the hipest ener^ channel. Above 1 Mev the integral 
flux for the hardest spectrum is of the order of 1,300/ (cm x see x steri 

For the other energy extreme we consider the flux at 50 Kev which gives 

6 5 2 

the values 5.5 x 10 and 7*8 x 10^/(cm*^ sec ster). At 400 Kev we have 

2 

4.5 and 4.6 x 10 , The data accumulation time was .00853 seconds as 

set hy telemetry requirements. We chose to set the "VTldths of the energy 
channels at approximately logarithmically increasing values to insure 
complete coverage of the energy range hy ovir nine charmel analyzer. The 
geometrical factor was designed basically to keep the largest expected 


TABIJS 2 , 1- Model Count Rates with F 


=^e*p{.E/E^) 


Sample energy 
interval (kev) 


Max. counts J =: 1. 6X10® J= 1.4X10® 1 = 1.0X10® J = 7. OXlo"^ . 

per channel = 50 kev EP*75kev Ei = 100 kev EP- 150 kev . 

o o ^ ^ 


40-60 

1023 

1011 

820 • 

518 

286 

60-90 

1023 

928 

885 

607 

362 

90-150 

1023 

789 

991 

783 

540 

150-220 

511 

256 

491 

479 

410 

220-340 

511 

76 

253 

330 

379 

340-600 

255 

8 

62 

131 

255 

600-1000 

255 

1 

2 

10 

51 

1000-1800 

255 

1 

1 

1 

4 

1800- CO 

255 

1 

1 

1 

1 
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flux in the lowest channel within the counting capability of the 
accumulators. The three lowest channels had 10 bit accumulators with a 
■mnvimuwi count rate then of 1023 counts per accumulation period. The 
next two accumulators had nine bits and the remaining four had eight 
bits. Table 2.1 demonstrates count rates for a set of hypothetical 
energy intervals and model spectra using a geometrical factor of .005 
cm ster. As can be seen, considerable latitude was available. This 
is particularly true for the higher energy channels in which we provided 
a much higher count rate capability than would be necessary for a pure 
exponential spectrum. 

It can be seen from equation 2.5 that with the geometrical factor 
chosen, there are still two free parameters,, which may be varied to 
determine the angular width of the collimator acceptance cone. It was 
anticipated that the rocket would have a roll rate of 5 roll/sec. or 
an angular rate of about 15 degrees per accumulation interval. 15 
degrees then became a reasonable angular width for the collimator. This 
width is indicated by the angle y in Figure 2.3. The remaining free 
parameter was adjusted to set the length of the collimator within the 
constraints of space in the rocket. Table 2.2 lists all the relevant 
dimensions of the collimator. 

Table 2.2 

Collimator Dimensions 
= .325 ± .0025 cm 
Rg = .ITO ± .0025 cm 
h = 2.1*67 ± .0127 cm 
Y - 15.0 ± .3 degrees 
G = .00l*8 ± .0003 cra^ster " " 


2. Deslgg of the Shielding 


There are several other Important design considerations 
for the collimator. Particles on the edge of the acceptance cone will 
scatter when hitting the collimator surfaces. These particles will have 
distorted energy spectra, as well as being outside the defined field of 
view. They will therefore introduce errors if they reach the detector 
and are counted. Particles stopping in the collimator wall will also 
produce Bremsstrahlung x-rays. Vnile the walls may be sufficiently 
thick to stop 3 or 4 Mev electrons they may not be thick enough to pre- 
vent significant x-ray count rates in the 3mm thick silicon detector. 

The scattering of electrons off the walls of the interior of the 
collimator cannnot be prevented. One attempts to reduce the problem 
by; (1) reducing the surfaces in direct line of sight of the detector 
surface and the entrance aperture, (2) providing scattering volinnes in 
which scattered particles may be reduced in energy and absorbed, and (3) 
orienting surfaces so that scattering particles are diverted toward the 
scattering volumes and away from the detector. In Figure 2.3 we show a 
diagram of the collimator in cross-section. The only surfaces in direct 
view of the detector surface and the entrance aperture are the sharp 
edges of the collimator baffles. These edges defined the acceptance 
cone and the geometrical factor. The acceptance cone was checked 
directly with energetic electrons. We will discuss this fully in 
section II. D. 

Again using the Bethe formula (eq. 2,3) we determined that 2.24 ram 
of brass was the average range of 3 Mev electrons. Adding 40% for 
straggling and insurance we decided that 3.2 mm of brass was sufficient 
for the outer layer of the colliinaton This was then backed with 3.9 mm 










2k 


of lead to attenuate the x~rays. It is itnportant to note that the 
bremsstrahlung production cross-section in lead Is about three times 
greater than in brass so that one prefers if possible to stop all the 
electrons in the brass. Direct line of sight between the lead and the 
detector was avoided to exclude x-rays produced by scattered electrons 
hitting the lead. 

!Ehe weight and size of the collimator were very important since it 
had to be carried on a rocket. We were forced therefore to keep the 
lead shielding to a bare minimum. We modeled the x-ray production in an 
attempt to determine this minimum shielding. We chose a very simple 
geometry for the task. A unidirectional flux of electrons of differ- 
ential spectra of the foinn of eq, 2.1 was assumed nomally incident on a 
plane target of brass and lead. A 3 mm thick silicon detector was 
assumed placed behind the target and the energy spectrum detected in the 
silicon was determined. This was then compared with the spectrum of the 
incident electrons. The calculation is described in detail in Appendix 
A-1. The results from the calculations for no lead are compared with 
the calculations for .39 cm lead in Figure 2.4. 

We have assumed that for the x-rays the detector had an effective 

2 

geometrical factor of Zircm star , where A is the area of the detector. 
It can be seen from Figure 2.4 that 3.9 mm of lead provided sufficient 
shielding for the case of a relatively hard spectrum. We used 3.9 mm of 
lead in the actual collimator and as will be seen later the spectrum 
encountered during the rocket flight was softer than this example. In 
section I. D we will compare these results with those obtained on the 
Goddard Space Flight Center Van de Graaff accelerator. " ~ 


RELATIVE DIFFERENTIAL SPECTRA 







The Aluminum Light Shield 


Since silicon detectors are highly light sensitive a light 
shield was included in the collimator to front the detector. We chose 
aluamlnum becaiJse of its high reflectance and opacity even into the ultra 



violet. Table 2.3 lists transmission coefficients for alimiinuin at 


several wavelengths. 

The aluminum foils for the windows were prepared by flash 
vacuum plating onto small glass slides that were dipped in a twenty per 
cent solution of polyvinyl pyrrolidone in ethyl alcohol. The thickness 
of the aluminum was determined with an interferometer. Each slide was 
examined for pin holes under an intense light. The regions with no pin 
holes were determined and etched around with a stylus. These regions 
were then floated off the slides by using distilled water to dissolve 
the polyvinyl pyrrolidone. The aperature of the collimator was placed 
under the foils and the distilled x^ater allowed to evaporate. The sur- 
faced adhesion of the aluminum to brass was sufficiently strong to with- 
stand the 34 g accelerations of the NlKe-Tomahawlc shake test, so no 
additional adhesive was required. 



TABLE .2, 3 - Transmission and. Reflectance of Aluminum Films at various 
wavelenghts, (taken from Hass and Waylonis, 1961) 


Wavelengths 


Film Thickness 

* 2200 A 

R% T% 

3000 A 
R% T% 

4000 A 
R% T% 

5460 A 
R% T% 

6500 A 
R% T% 

3200 A 

88.5 

3.2 

90.0 

1.8 

91. 1 

1. 1 

90.4 

0,5 

89.6 

0.4 

4000 A 

90.6 

1.1 

91.4 

0. 5 

92. 1 

0.4 

91.2 

0.2 

90.3 

0.2 

5000 A 

91. 5 

0.3 

92.0 

0. 1 • 

92.5 

<.l 

91. 5 

<. 1 

90.6 

<1 


Iff »Wn 


28 


II, C. Effects Of The Light Shield On Bnergy Resolution 

2 

10,000 ^ of aluminum is a dead layer of 270 ygm/cm « Figure 
2,5 gives the energy loss of various energy electrons in this thickness 
of aluminum. A 30 Kev electron loses about 2 Kev. The energy loss be- 
comes even less as the energy of the electrons is Increased, so that for 
the largest portion of the energy region of interest, the foil causes 
little distortion in the average energy. There are other possible dis- 
tortions that should be considered as a slight digression into the energy 
loss process lirLll revival. 

As an electron traverses a material like aluminum it loses energy 
by numerous ionizing (i.e. inelastic) collisions with atomic electrons. 
Bie vast majority of these involve almost no angular deflections (Bethe 
1936, Landau 1946). In addition the electron undergoes elastic col- 
lisions with the screened nuclei of atoms in which it may suffer sizeable 
angular deflections. Radiative energy loss due to production of 
brerosstrahlung is unimportant in both collisional processes for primary 
electron energies below 10 Mev. (At a primary electron energy of 1.5 
Mev the ratio of collisional to radiative energy loss in aluminum is 
6.9 X 10 (Berger and Seltzer, 1964)X Since the energy loss and the 
angular scattering are essentially stochastic processes, the electrons 
emerge from a thin foil with a distribution of energies and directions. 
Hie propensity of electrons to lose a varying amount of energy over any 
short distance is usually referred to as energy straggling and it is well 
described by the Landau distribution. (Deamaley and Northrop, 1966). 

We see then that unidirectional monoenergetic incident particles emerge 
from a thin foil with a distribution in both energy and direction. Some 
particles are scattered back in the initial direction and lost completely. 
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Figure 2.5. Energy loss of electrons in 270 usn/cm of aluminum 
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We are therefore required to examine the transmission, probability as 
well as angvilar and energy distributions. 

Experimental and theoretical data does exist over the energy range 
of interest although the coverage is not very complete. The data that 
does exist indicates that extrapolation to energies not specifically 
covered can be done with a fairly high degree of confidence. 

Ranter (1961) has determined the transmission coefficient, Ti^ , 
for electrons in the energy range 1 to 10 Rev through various thickness 
of aluminum foil. The beams as a function of reduced energy, E , were 
normally incident and the transmitted particles were measured over 2fr 
steradians. The transmission coefficient may be expressed as 

a transmitted particles over 2ir steradian (2.T) 

^T incident number of particles 

The coefficient is normally scaled in terms of the reduced energy 
E which is given by the relation E = E/E^ with E as the energy of 
the Incident electrons and as tlie penetration energy of the foil as 

defined in appendix A-2. Figure 2.6a is a slightly modified version of 
Ranter's results for aluminum. It has been modified for easier com- 
parison with Figure 2.6b which is for the energy range 500 Rev to 5 Mev 
determined from the results of Berger and Seltzer (1968). The agreement 
is quite good even though the energy ranges are very different. In 
Figure 2.6c is a curve determined from the empirical formula of Robetich 
and Rats (1968), The empirical formula is very simple to use since it 
can yield values at all energy without scaling from a graph. For this 
reason we used Robetich and Rats's results. For AO Rev electrons, E 

O 

is about A for a 10,000 A aluminum foil. It is evident that from 


FRACTION TRANSMITTED 












Figure 2«6c that the correction is quite insignificant over most of the 
energy range of interest. 

Figure 2.7 shows energy spectra of transmitted particles for 
incident beams of monoenergetic particles scaled at two values of E 
again taken from Berger and Seltzer (1968). Interpolating we can de- 
termine that for an incident energy of 40 Kev the resolution (Full-Width- 
Half-Max) is of the order of 3.5 Kev. This resolution obviously improves 
with increasing energy. As we will see this is small compared to the 
detector-amplifier energy width. 

The final problem associated with the foil is the angular scatter- 
ing and like the other problems is most pronounced at lower energy. 

Figure 2.8 shows a plot from Kanter (1961) of the most probably scatter- 
ing angle as a function of the reduced energy, E . For E = 4 (about 
40 ICev) the most probably scattering angle is less than fifteen degrees 
from the normal, but some 5% of the particles have scattering angles 
greater than 50 degrees. As the solid state detectors have a signifi*^ 
cantly larger dead- layer off to one side (see Fig. 2.1) something on 
the order of 2% of the lowest energy counts may be lost. 




37 


II, D, Modeling The Spectrometer Response Functions 

For energies below 10 Mev electrons lose energy mainly in ioniz- 
ing collisions with atomic electrons. Secondary and tertiary particles 
are produced with most having kinetic energies smaller than the average 
ionization potential of the absorber atoms. (Deamaley and Northrop, 
1966) . The cascading particles continue the process of impact ionization 

until their kinetic energies are all below E , the energy necessary to 

S 

raise one electron out of the valence band of the silicon into the 
conduction band. (Taylor, 1963). There are several processes other than 
production of conduction electrons and holes which act as final reposi- 
tories for energy so that on the average the final energy to charge pair 

ratio is slightly higher than E . E for silicon is 1.08 ev per 

S S 

electron-hole pair (Deamaley and Northrop, 1966) while the measured 
mean energy per electron-hole pair produced by ionizing radiation is 
W = 3.6 ev. Neglecting any charge loss processes we would expect 

Q = Ee/W (2.8) 

to be collected on the electrodes. E is the energy of the incident 
particle and e is the electron charge. For a 100 Kev electron stop- 
ping in silicon this formula yields 5.3 x 10 coulombs. This means 
we require a very sensitive charge integrating preamplifier as the first 
stage of £in amplification system. 

As we have indicated the average energy loss per path length for 
electrons stopping in different materials varies inversely with energy 
below an Mev. This means that energy deposition or ionization is an 
Increasing function as an electron penetrates into a silicon detector. 







A significant quantity of collectable charge is thereby produced very 
near to the end of the particle path (Kobetich and Katz, 1968) . The 
collection time of the charge or equivalently the pulse rise time of 
detector signal is a function of the particle energy. With this fact 
In mind one normally designs pulse circuits with time constants longer 
than the longest expected pulse rise times. This provides for complete 
charge collection and a linear relationship between charge and particle 
energy. If one chooses to make the pulse circuit time constants of the 
same order as the charge collection times, thereby gaining increased 
count-rate ability, the collected charge and particle energy are no 
longer linearly related and the energy resolution is degraded at high 
energies. With the wide d 3 mamic range of predicted particle fluxes in 
REP events, our main concern was for a great count rate capability. As 
we also did not Intend to take advantage of the fine energy resolution 
at high energy available from the P-I-N detectors, short pxxlse circuit 
time constants seemed to provide a very reasonable compromise for our 
experiment. 

The charge sensitive preamplifier and post- amplifier combination 
were to have a maximum average count rate of 200 Khz. Allowing for 
random spacing between counts we set a minimum time resolution of the 
order 500 nsec. In Figure 2., 9 we shot^ a graph of positive charge 
collection time versus particle energy for our three mm thick detectors. 
For energies in excess of 700 Kev the collection times are as long or 
longer than the 500 nsec time constants, thus giving us a non-linear 
pulse height vs. energy relation. 

In Figure 2,10 we show the basic amplifier system.' The Initial 
stage is a charge sensitive preamplifier, which is followed by two post 
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Figure 2 , 10 , The basic ahpilfier system 
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amplifiers in parallel with each other. The pulse height discriminators 
are shown in Figure 2.11, The five lowest discriminators received the 
pulse signal from postamplifier A and the four highest received the 
pulse signal from postamplifier B. The discriminator logic allowed only 
one counter to receive a pulse for each charge pulse on the detector. 

The discriminators were also Inhibited during the time a pulse was being 
processed. The charge terminator at the input of the preamplifier en- 
abled us to calibrate a fast rise time (less than 10 nsec.) pulse gen- 
erator in terms of equivalent energy pulse. As the rise time of these 
pulses were very short compared to 500 nsec, we have, therefore, a 
source of equivalent energy pulses that would yield a linear output pulse 
height versus equivalent energy (charge) input. 

Also shown in Figure 2,11 are the detector bias current and 
temperature monitor. The normal detector reverse bias currents were of 
the order of a few microamps. Any large deviations were indicative of 
break-do^m. The over-all noise of the detector amplifier system was 
quite temperature dependent so that significant degradation of resolu- 
tion, would occur at high temperatures. The monitors were included as 
warnings of these possible malfunctions. 

As we have indicated, the output pulse height of the amplifier 
systems as a function of input pulse height and rise times is a very 
important characteristic of this system. Figure 2,12 shows a typical 
calibration determined with a variable rise time linear pulser. t'Jhat 
is displayed is the amplifier efficiency, H , defined by the relation 

gain at rise time t _ _ , 

gain at rise time of 10 nsec 


E 
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The amplifiers were perfectly linear as a function of input pulse- 
height so that these curves were independent of pulse-height. As we 
see the efficiency begins to fall off at about 100 nsec rise time, 
which is one-fifth of the preamplifier time constant. 

The resistor string across the top of Figure 2.11 is a simple 
voltage divider. The relative values of these voltages determine the 
level at which each discriminator triggers. If the system were linear 
it would be necessary only to make these voltages proportional to the 
desired energy thresholds. Since the output pulse height was not 
linearly related to the energy of the Incident electrons we had to use 
a somewhat more involved procedure: 

For eacli preamplifier-post amplifier combination we determined the 
efficiency as a function of rise time (Figure 2,12), Knowing the charge 
collection time as a function of energy we then produced the output pulse 
height as a fimction of incident energy by performing the convolution of 
E with the collection time. The collection time was determined by cal- 
culating the time reqtiired for the most distant charges produced by the 
ionizing electrons to be swept to the electrodes of the detector. The 
entrance surface for the electrons was the p-type gold surface of the 
detector. As the junction is maintained at a reverse bias the positive 
carriers (holes) are collected at the gold electrode and the negative 
carriers (electrons) are collected at the aluminum electrode. Ihe total 
rise time of the electron charge collected on the aluminum electrode is 
is therefore independent of energy since total collection requires the 
electrons be swepted across the entire thickness of the detector. The 
rise time of the hole charge on the gold electrode Is given approximately 
by the expression 
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Figure 2.12. Amplifier output pulse height at constant input pulse height as a function of 
input pulse rise time. 
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vhere d is the detector thickness, r the depth of penccratlon of the 
ionizing particle, is the hole mobility and V Is the potential 

between the electrodes (Deamaley and Northrup, 1966), As shown by 
Weber (1964) the average penetration depth of electrons into various 
materials is well represented by the formula 


where E is the energy of the incident electron and A, B, and C are 
constants characteristic of the stopping medium. This formula was 
applied to the graph of penetration depth vs. energy given by Deamaley 
and Northrup (Deamaley and Northrup, 1966) with the resulting values of 


A = 2.205 X 10 
B = .9476 
C = 2,292 X lO' 



Combining equations 2.9 and 2.10 we produced the curve of hole charge 
collection time vs. energy as shown in Figure 2.9. This result was 
then folded through the efficiency curves for each preamplifier- 
postamplifier combination to produce the pulse-height vs. energy curve 
shown in Figure 2.13, The pulse height is expressed in equivalent 
energy units. The relationship is linear up to an electron energy of 
320 Kev at which energy the hole rise time becomes greater than 0,1 psec 
and part of the hole charge is lost. The vertical lines superimposed on 


PULSE HEIGHT (KEV) 



ELECTRON ENERGY (KEV) 

Figure 2.13> Pulse height as a function of energy. 

The solid line Is the theoretically determined pulse height. The vertlcle lines were determined 
experlmentolly. a* 



the curve are the pulse helots measured from electrons produced on the 
2 Mev Van de Graaff accelerator at Goddard Space Flight Center. We will 
discuss this further in the next section. 

Ihe values of the resistors in the pulse height discriminator were 
now determined (in a straight forward manner) from this curve. The 
final equivalent energy triggering level of each discriminator was 
checked with the 10 nsec rise time pulses from the linear pulser. 
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II, E. Calibrations 

Ibe calibration of our detectors and associated electronics was 
performed on the 2 Mev Van de Graaff and 100 Kev electron-proton accel- 
erators at Goddard Space Flight Center- These calibrations were carried 
out in two stages. The first stage was a determination of the optimum 
performance characteristics and dead- layers of our lithium drifted solid 
state detectors using electronics with relatively long circuit time con- 
stants and electron beams of low flux. The pulse height analysis was 
performed by a 10,000 channel Nuclear Data Pulse Height analyzer. The 
second stage was the calibration of the completed flight model nine 
channel electron energy spectrometer with collimator and altaninum window. 

1. The Determination Of Optimum Performance 

The detectors were placed in the scattering chamber on a 

moveable table which cwsuld be positioned by remote control to allow for 

2 

small rotations and lateral movements. A 5 cm lithium drifted solid 
state detector which could be moved in and out of the electron beam was 
positioned between our detectors and the chamber entrance. We used this 
••monitor" as a check on the beam intensity. The electrons were magnet- 
ically analyzed before injection into the scattering chamber and there- 
fore had energy resolution far better than the 25 Kev resolution of the 
system being calibrated (on the order of ,1%). 

The pulse height spectrum for each sampled detector was taken 
at numerous energies betoeen 100 Kev and 2,0 Hev, A typical spectrum is 
sho\m in Figure 2.14a for a 150 Kev bcaia. The large gausslan shapv'd 
peak is called the absorption peak and is the result of particles leav- 
ing all their energy in the active region of the detector. The long 



tail extending toward lower energies results mainly from particles escap- 
ing out the front of the detector and a small portion out through the 
sides. This is termed the escape component. The increase toward lower 
energy (near 50 Kev) is due to noise counts above the lower threshold of 
the analyzer. 

The monitor detector, like the ones being calibrated, has a 
lower threshc d and a thin gold layer on the entrance surface of about 
20 Uffa/cm . This meant that ve could not rely on the monitor to obtain 
absolute covinting efficiency. Though the corrections were generally 
small (less than 105?) we ad;Justed efficiencies to correspond with the 
tabulations in Berger and Seltzer (1969). They have written a monte 
carlo coB^Juter program which follows the lif» histories of 10,000 
electrons and their secondaries as they stop in ve.rious materials. 

Berger and Seltzer (1969) contains complete results of energy deposition 
spectra for electrons stopping in silicon and comparisons of the 
theoretical results with numerous experiments for the energy range 150 
Kev to 10 Mev. For each energy of interest we produced the theoretical 
spectra according to the procediu’es outlined in Berger and Seltzer 
(1969). Such a spectrum is shown in Figure 2.15b. V?e then normalized the 
measured spectra by setting the area within one full-width-half-maximum 
of the absorption peak equal to the same area tinder the theoretically 
determined spectra. The normalized spectra corresponding to that in 
Figure 2.lU is shown in Figure 2.16. For energies greater than 2 Mev 
for which we could not obtain actual electron measurements we adopted 
the spectra from Berger and Seltzer (1969). 

The spectra determined on the Van de Graaff, renormalized as we 
have explained and corrected for transmission through the aluminum foil 
gave us a set of functions F(e,E) , where e is the abscissa in 


XlO^ELECTRONS/KEV 
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Figure 2,15 and E is the energy of the incident electrons. To de- 
termine the response of a channel from the nine channel analyzer with 
Icwer threshold and upper threshold e^ to electrons of energy E 

we evaluate the integral 



F(e,E) de 


( 2 . 11 ) 


The nine response functions for* two spectrometers are shown in Pigs. 
2.36a & h by the solid lines. We built six complete systems and deter- 
mined the response functions in the above manner for each system. 

2. Final Calibrations 


During the second stage of calibrations there were five specific 
characteristics of the coB^jleted systems that we checked: (l) the non- 

linear pulse hei^t vs. energy relation discussed in section II. D, (2) 
the energy response functions of the nine channel spectrometers, (3) the 
dependence of the response functions on particle coui*t rates, (U) the 
efficiency with which x-ray counts were excluded by the collimator 
shielding, and (5) angular acceptance of the collimator. 

The collimator-preamplifier assembly with the aluminum light 

2 

shield was mounted on a remotely controlled moveable table with a 5 cm 
solid state detector again used as a ’’monitor". To reduce scattering in 
the collimator the electron beam entered the scattering chamher through 
a 2 mm diameter aperture and the table position was adjusted until a 
maximum counting rate was achieved. 

The pulse height vs, energy relationship was determined by using 
electrons of eighteen different energies between 50 Kev and 2 Mev, The 


i 

li 


i! 

I 
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output from the postamplifiers was displayed on an oscilloscope while 
electrons were incident on the detector and a calibrated linear pulser 
provided a reference signal to determine the equivalent energy pulse 
height. The results for one such calibration are shown by the crossed 
lines in Figure 2.13. The agreement between the ’modeled relation and 
the measured one were in each case of the same quality as shown in this 
example . 

To determine the response functions of the nine channel analyzer 
required IjO energies between 20 Kev and 2 Mev, The beam intensity was 
checKed with the monitor detector before and after each energy. If the 
intensity of the beam varied by more than ten per cent the data were re~ 
taken. To deteamiine the absolute efficiency the count rate in the 

2 

monitor detector ■vras corrected for transmission throu^ 20 pgm/cm of 
gold (the entrance electrode of the monitor detector) and for counts be- 
low the noise threshold (tising the data from Berger and Seltzer (1969)). 
The crossed lines in Figure 2.16 are the points determined in this 
manner. The error indicated by these points is due strictly to the 
variations in the beam intensity over the comting period since each 
point was determined by well over 20,000 counts. The response function 
as determined by the integrals in equation 2.11 were in excellent agree- 
ment with the points determined in this calibration procedure, as is 
evident from Figure 2.16. 

For selected energies between 50 Kev and 2 Mev, points on the 
response functions were also checked at six different beam intensities. 
The intensities were adjusted so as to give approximate count rates of 
10, 25, 50, 100, 150 and 200 kilohertz in the detectors". "The values for 
several energies and different coimt rates for one sucli calibration are 
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shown in table 2.4. In each case there was no significant change in the 
value of the points determined. This res-ult implied that the pulse 
circtdtiy introduced no count rate dependence and that there was no 
pulse pile-up at even the high rate of 200 kilohertz. It was to this 
end that we designed the pulse circuitry with such short time constants. 

The collimator opacity to x-rays was the next characteristic 
that was checked. The first element of the collimator baffles was re- 
moved and replaced by an identical piece with no entrance aperture. The 
coUimiiitor was then subjected to beams of electrons of varying intensity 
and several energies between 50 Kev and 2 Mev. Table 2.5 shows the re- 
sults of this determination. In each case the x-ray count rates are 
very small so that we could ignore them when performing the final 
analysis of the data. 

The geometrical factor 6 for the collimator was determined hy 
the parameters , Rg and h as defined in figure 2.2 and related to 
G by equation 2.5. As it is very difficult to produce an isotropic 
hesim of energetic electrons we had to attempt to determine the parameters 
by a less direct approach then measuring G itself. We constructed a 
mounting which allowed us to rotate the collimator-detector housing 
about an axis that passed directly throxigh the center of the exposed 
face of the solid state detector. The detector and collimator were 
e^qposed to a monodirectional heam of electrons while supported on the 
rotatable table. The beam was diffuse and covered the entire aperture 
of the collimator This assembly was then rotated in steps of one de- 
gree about the axis while the count rate was monitored. In this fashion 
we determined the angular response of the collimator to a monodirectional 
beam making on angle a with the collimator axis. This response is 
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TABLE 2. 4 - Uncorrected counting efficiencies at several rates for 
discriminator 3. (GSFC - April 14, 1972) 


Approximate particle rate at detector 


E=100 kev 

lOHz 

50Hz 

150Hz 

200Hz 

CH 1 

. 115 

. 116 

. 115 

. 115 

CH - 2 

.775 

.775 

.773 

.773 

CH - 3 

.021 

,022 

.023 

.021 

CH - 4 

0.0 

0.0 

0.0 

0. 0 

CH - 5 

0.0 

0,0 

0.0 

0.0 

CH - 6 

0. 0 

0. 0 

0.0 

0.0 

CH - 7 

0. 0 

0.0 

0.0 

0, 0 

CH - 8 

0. 0 

0. 0 

0.0 

0. 0 

CH - 9 

0.0 

0.0 

0. 0 

0. 0 


E=300 kev 

lOHz 

SOHz 

150HZ 

200Hz 

CH - 1 

.0179 

.0173 

.0180 

. 0175 

CH ~ 2 

.0243 

.0242 

. 0242 

.0241 

CH - 3 

.0517 

.0517 

. 0515 

. 0520 

CH - 4 

.735 

.736 

.736 

.738 

CH - 5 

. 152 

. 155 

. 154 

. 163 

CH - 6 

0. 0 

0.0 

0.0 

0. 0 

CH - 7 

0.0 

0.0 

0.0 

0.0 

CH - 8 

0.0 

0. 0 

0.0 

0.0 

CH - 9 

0.0 

0.0 

0.0 

0.0 

E-1200 kev 

lOHz 

SOHz 

150Hz 

200Hz 

CH - 1 

2, OOE-4 

2. 07E-4 

2. 03E-4 

2. 03E-4 

CH - 2 

4. 27E-4 

4. 27E-4 

4. 33E-4 

4. 26E-4 

CH - 3 

2. 19E-3 

2. 18E-3 

2. 17E-3 

2. 21E-3 

CH - 4 

6. 58E-3 

6. 57E-3 

6. 62E-3 

6. 57E-3 

CH - 5 

3. 95E-3 

3.99E-3 

3. 94E-3 

3, 98E-3 

CH - 6 

6. 8 IE -2 

6. 81E-2 

6. 78E-2 

6. 79E-2 

CH - 7 

1. 95E-1 

1.96E-1 

1. 91E-1 

1. 92E-1 

CH - 8 

7. 59E-1 

7. 51E-1 

7. 62E-1 

7. 59E-1 

CH - 9 

Q. 0 

0.0 

0. 0 

0.0 



TABLE 2, 5 - X-ray c< unts per 


Energy (kcv) 

CH-1 

CH-2 

CH-3 

CH-4 

CH-5 

400 

p. 0 

0. 0 

0. 0 

0. 0 

.20E-6 

600 

0.0 

0. 0 

0, 0 

.43E-6 .62E-6 

800 

0. 0 

0. 0 

0. 0 

.71E-6 .llE-5 

1000 

0. 0 

0, 0 

0. 0 

. 12E-5 . 14E-5 

1500 

0.0 ■ 

0, 0 

. 19E-6 

.30E-5 .35E-5 


1 

i 




CH-6 

CH-7 

CH-8 

CH-9 

.21E-4 

0. 0 

0.0 

0.0 

.31E-5 .40E-6 

0.0 

0.0 

.69E-5 .41E-5 

0.0 

0.0 

.94E-5 .77E-5 

.29E-6 

o 

o 

.27E-4 .30E-4 

. lOE-4 

0. 0 



6o 


shewn in Figure 2,17 by the short horizontal lines. Using the values of 
, Rg and h given in table 2.2 ve deterEiined numerically what this 
angular response should be. (This calculation is described in appendix 
A- 3). These values are shown by the solid line in Figure 2.17. As can 
be seen from the figure the measured width of the collimator is slightly 
broader than the width determined from the measured values of the para- 
meters but well within the error bars. This excess rai^t be due to a 
small amount of elastic scattering from the sharp edges of the baffles. 
Increasing the values of or Rg , or decreasing the value of h will 

not account for this slight bowing, in that the meastired angular width 
goes to zero at the predicted value of a . 


RESPONSE 
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II. F. Energy Response Of The Analyzers To Protons 

The ptO-se height coming from the peramplifier and postaraplifier 
combination, as we have indicated, was a function of the particle energy. 
This was due to the dependence of the charge collection time on the 
depth of penetration of the ionizing paidiicles into the silicon de- 
tector. The particle penetration, depth is a function of the particle's 
mass and charge as well as energy. One should therefore anticipate that 
the eq,xiivalent energy thresholds of the analyzer would he charge species 
dependent. Though it was anticipated that the primary somrce of ioni- 
zation would he electrons, we had to he prepared to identify protons if 
they were present and compensate for them in the data analysis. There 
was a separate proton detector included in the rocket instrumentation. 

Protons straggle far less than electrons because of their greater 
mass. In addition a one Mev proton ifill penetrate only on the order of 
17 microns into silicon. (Northcliffe and Schilling, 1970). At a depth 
of 17 microns the hole charge was completely collected so that the re- 
sponse of the energy analyzer was linear to the proton energies of in- 

2 

terest. The effect of the 270 ygm/cm of aluminum in the light shield 
2 

and the 31 ygm/cm of gold on the entrance electrode was deteimned by 
the use of the stopping power tables in Northcliffe and Schilling (1970). 
As in the case of the determination of the electron response functions 
we took into account the charge collection deficit of the negative 
carrier while assximing the holes were completely collected (see section 
2-D), Table 2.6 lists the equivalent energy thresholds of one of the 
electron spectrometers. .. - 


TABLE 2. 6 - Energy Response of 45~degree spectrometer 
to Protons, 


Proton Energy (kev) 


CH-1 

173 - 185 

CH-2 

185 - 214 

CH-3 

214 - 270 

CH-4 

270 - 339 

CH-5 

339 - 364 

CH-6 ■ 

364 - 506 

CH-7 

506 - 736 

CH-8 

736 - 1320 


CH-9 


1320 - oo 
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II ’G* PCM Telemetry 

The data from the rocket was retviriied in Pulse Code Modulated (PCM) 
binary fonaat. The carrier frequency was 2li0.2 MHs and the bit rate 

C 

was 1.2 X lO"^ bits/sec. The bit stream was subdivided into 8-bit words, 
32-word frames and 32-frarae master frames. A single 8-bit word was 
66.1* ]isec long. All three of the spectrometers and the BTL proton 
experiment accumulated data simultaneously over an accumulation time of 
8.53 msec. At the end of this period the counts were stored in a buffer 
and read serially into the bit stream while the spectrometers were 
accumulating fresh data. This accumulation period was four frames (32 
words) long- The spectrometers each required 10 words while the BTL 
proton detector required nine words. The data from each of these 
detectors was therefore injected into the bit stream once every four 
frames using the same 10 data words each frame as sho*<jn in table 2.7* 
This is referred to as subcommutating. The tenth word on the BTL proton 
experiment was set at 11111111 to serve as a synch-word and thereby 
identify the position in the four detector cycle. The analog monitoring 
data (usually referred to as housekeeping data) such as the detector 
temperature and current monitors was converted to digital (A to D) 
and suhcommutated once every 32 frames into the bit stream. Table 2.7 
lists the basic telemetry frame. 


TABLE E. 7-Telemetry Channel Assignment F-34 
Video: Frequency E40.2 Mhz 

Information rate: 8 bits/word X 32 words/frameX 469 frames/ 

is 120, 064 bit/sec 

Each word 66. 4 microseconds 
Each bit 8, 3 microseconds 


1 I I 0 0 0 10 synch-word 

0 1 0 0 0 0 1 I synch-word 

Q 0 10 o' I 0 0 0 0 open 

sub-corn digital 2 

sub-corn analog data 

Faraday receivers 15.011 MHz 
Search coil 

Faraday receivers 3. 883 MHz 
Channeltrons protons 0 degrees 
Channeltrons electrons 45 degrees 
Channeltrons electrons 90 degrees 
UOM and BTL electron- oroton experiments 


Search coil 


Gerdien orobe 


NDRE Sun sensor 


NDRE Lvman-Alpha probe 








CHAPTER III 


THE PLICfHT PARAtffiTERS OF ROCKET F-3*f 
III. A. The Tra.jectory 

The rocke'l?s trajectory was determined by a slant range system. The 
actual trajectoiy parameters were prepared by the ballistics section of 
the Andoya Rocket Range, NDRE. In Figure 3.1 we show a diagram of this 
trajectory. The rocket attained a maximnn altitude of 15^.43 km at 197 
seconds after lift-off. The average flight azimuth was 320 degrees. The 
rocket attained 90 km altit\ide at 80 seconds into the flight and remained 
above this altitude until 314 seconds. For the portion of the flight 
above 90 km the rocket traveled 113.5 i .9 km horizontally at azimuth 
320 degrees which is some 86 km north and 73 km west. 

The aspect determination depended upon the direction of the mag- 
netic field and the direction of the sun. The magnetic field direction 
was taken from Wiley and Barish (1970 ) . Prom this model it was determined 
that the azimuth of the magnetic field varied from 0 degrees to - 1.3 
degrees over the entire trajectory, while the dip angle varied from 
72.5® to 72.8°, The solar elevation was 28.9® and the solar azimuth was 
102.4° at the time of launch. The solar position would have changed by 
less than .3° over 400 seconds. We therefore took the magnetic field 
direction and the sun direction as constants for the aspect determination. 

We also note at this point that according to Stassinopoulos (1970) 
the magnetic field line at 100 km at the position of the rocket extended 
out to 6.4 Earth radii at the equatorial plane. - - 
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III.B. Determining Rocket Aspect 


The flux of the precipitating particles must he determined with re- 
spect to a reference system fixed on the Earth. The obvious choice is 
one with the Z-axis along the local magnetic field of the Eaiiih since 
many of the properties we wish to deteimine relate directly to the mag- 
netic field. The rocket carried an aspect magnetometer, which determined 
the angle between the rocket axis and the magnetic field as well as giv- 
ing an accurate account of the rockeifs spin rate. In addition, there wa.s 
a solar aspect sensor which determined the angle between the rocket axis 
and the sun direction as well as giving a measure of the spin rate. 

Knowing the elevation and azimuth of the sun and the direction of the 
local magnetic field it is possible with the data from the sxjn sensor 
and the magnetometer to determine the orientation of the rocket in space 
at all times. 

Unfortunately after JO seconds into the flight the magnetometer data 
were lost. After a number of unsuccessful attempts we discovered that 
the electron data from the downward looking spectrometer could be "used to 
determine the information lost by the magnetometer. The procedure was 
somewhat involved but many of the conclusions of this research are based 
upon the reliability of the procedure so it warrants detailed discussion. 

Figure 3.2 shows the angle between the noarmal to the rockeu axis 
and the sun direction as a function of time determined from the sun 
sensor data. (E. Thrane, private communication). The altitude of the 
rocket from 85 seconds to 300 seconds is above 90 kilometers where the pres- 
sure of less than .75 microns is small enough that the drag due to air be 
neglected and we may assiime that the rocket was in free fall. The nose cone 
was spin-balaiiced about the rocket axis so that this axis was a rotational 





syiranetry axis. The free fall motion then consisted of a rotation about 
the symmetry axis at on angular frequency w and a precession of the spin 
axis about the direction of the angular momentun at an angular frequency 
Si. (Goldstein, 1950). 

To determine the eqtiation of motion of the rocket it is con- 
ceptually simplest to use three coordinate systems which have the same 
origin and are related by three sets of Eviler angles. 

1. The rocket coordinate system with 7 .^ along the 

rocket axis, along the magnetometer axis and 

Y chosen to form a ri^t handed set. 

2. The angular momentum system with along the 

direction of the total angular momentum and 
and Y fixed within a right hand system as 

Jj 

explained below. 

3. The sun-rocket system with Zg along the sun 
direction, 3^ in the local Earth horizontal plane 
and Y chosen again to form a ri^t handed set. 

We are interested here in only the 3 degrees of rotational freedom 
and will therefore have to fix some of the additional coordinates we are 
creating. 

Figure 3.3 shows the choice of Euler angles that transform from 
coordinates 1 to coordinate 2. The transfori!*ations consist of three 


simple rotations. One in the ^ " plane, followed by a rotation in 

t 

the Y - Z.^ plane and the last in the X* - Y” plane. Angle a is 

ll 

t 

the half coning angle of the rocket precession about Z . X labels 

h 

the line of nodes which moves aroimd the origin in the - Yj^ plane 
as the rocket processes so that angle ^ is the angle of precession. 


Angle is the rocket roll angle and is measured from the moving line 


of nodes. 
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In Figure 3.U we show the choice of EuXer angles that transform 
from • coordinates 2 to coordinates 3. We have set the angle between 
X-. and the line of nodes to zero removing one of the excess degrees of 

ij 

freedom. This fixes the direction of 35^ and , 

The unit vector pointing along the rocket axis in the rocket system 
is 



Bf performing the successive transformation we express this in the sun- 
rocket system 3 


rcosijj sina + sint{)(cosy cos C sina+sinj^ cosa) 
ft = I sin^ sin C sina - cosij>(cosu cos 5 sina+siny 
i - siny cos C sina + onsn . 


sina+siny cosa) 


(3.1) 


The sunsensor look direction, the rocket ••’xis end the sun direction 
are related as shown in Figure 3.5 

The angle X » which is the angle between the sun direction Zg and the 
normal to the rocket axis (sunsensor look direction), is measured once 
per roll by the sunsensor. The angle X between the rocket axis and 
Zg is related to the projection of Jig on Zg by the expression 


')■ 

j 
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sinX = a cos 5 + b *= cosX 


(3.2) 


with 


a -■ - sinp sina 
b = cosli cosa 


T 

(3.3) 


We are interested in determining the angles V and a which may be 
expressed in terms of a and b by the relations 



(3.U) 


We win take 
quency and ^ is 
nodes. The period 
from which in turn 
a and b by 


5 - f2t + ^ where is the precession angular fre- 
an initial angle between the axis and •*'he line of 

T = 2TT/fi is determined from the data in Figure 3.2 
we determine . The values of sirQ( are related to 


sinY = a + b 
'^max 

- a + b 


(3.5) 


To take advantage of the large number of data points determined from 
the sunsensor, rather than use relations (3. 5)» we least-squares fit the 
equation (3.3) to the data in Figure 3.2, solving for the constants a 
and b. Equations (3.4) then determined the two possible sets of solutions 
for a and b sho\m in table 2.1. The second set of solutions would imply 
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that the rocket turned nose dovm for part of the precession and nose up 
for part of the precession. The transversal magnetometer continued to 
function until shorbly before the doors of the rocket were opened. The 
indications were that the rocket had only a small precession cone before 
the magnetometer began to misf unction at 57*1 seconds. The doors opened 
at 57*6 seconds which is the time when the solar aspect sensor began to 
function. It would have required a sizeable impiolse to change the rock- 
ets orientation from the small precession cone at 57*1 seconds to the 
very large cone indicated by solution set 2 in table 3.1. As we will 
discuss in more detail later» the second solution may also be partly ex- 
cluded by an examination of the particle data. 

Table 3.1 

1 2 
]X = 78.21® +2® p = 173.96® + 2® 

a = 6.0U® + .5® a = 101.79° + .5® 

Of the initial five unknown parameters ? , 4^ , a , p and l|> we 
have determined three: C > w and p . <{i is the rocket roll angle with 

respect to the line of nodes . We take <p = on, "*■ and determine u 
by the following considerations. The sunsensor recorded the time between 
the crossings of the simsensor look direction and the sun direction. As 
the line of nodes Eiakes one 36O® precession every precession period the 
simsensor data recoided one more sun crossing than the number of roils - 
We simply coiint the nxmiber of crossing over the total of full precessions 
and subtract the number of precessions. This result is divided into the 
time to determine the roll period v^’ich is then converted to a frequency. 
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The value of u determined by this procedure was identical to that de- 
termined from the auto-correlation analysis of the particle data which 
will be explained in section 3.E. 

The determination of was somewhat more arduoias and proceeded 
as follows. The staisensor look direction is given in the rocket co^ 
ordinate system 1 by 



In the sm rocket system (3) this becomes 



”cosi|>(cos§cos4>-sin 5:osO!sin4>) - sinij?{cosu(sin§ 20 s({i+cos §:osasin<})) 

- sinvisinasin4>} 

sini{/(cos§cos4i-sin5:osasin4i) + cosi{»{cosii(sin5:os4i+co3 5cosasin4)) 

- sinjisinasintj)} 

sin}i{sin^os^cos 52 osasin<f») -j- cosMsinasin<Ji 


A 

The Z-component of Ugg was a maximum vxhen the sun direction crossed 
the look direction of the sunsensor. The sunsensor sees the sun when 



dt 


0 


Solving for we find 


(3.6) 


<j) = - cot + tan 

Tq ^ 




*tnnWcos(i}t+f; )(n+wcosa) + cjsina 

o 


tan|Jsin(J2t+ t.Q){o)+52cor.a) 




(3.7) 
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The numerous values of t were then substituted in to equation 
( 3 . 7 ) determining a set of values of 4*^ which were in turn averaged 
to determine the final value of 4*^ • As it happened no single value of 
4>^ determined by equation (3.7) differed from 6 q more than 

.005®. 

Figure 3.6 shows to what extent the rocket aspect is determined by 

the angles so far determined. The rocket axis moves on the surface of a 

cone. The axis of the cone is titled at an angle a with respect to 

the Z„~axis. The angle 4^ » which is the last unknown angle, is the 
b 

azimuth of the cone in the Xg - Yg plane. 

As the angle a fixes the angle that L makes with the X_ - Y_ 

O D 

plane, we see from Figiire 3.6 that the angle 6 between L and B is 
fixed once 4 ^ is determined. The angle 6 was to be determined 
directly from the magnetometer data. Since t’ > magnetometer failed we 
resorted to using the particle data itself to fix the value of ip which 
in turn fixed 6 , 

To elucidate, in Figure 4,1 we show the count rates of several 
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channels of the i*5°» 90® and 135® detectors. The 135® detector was to 
look.jnainly down the field line, the 90® detector was to look mainly per- 
pendicular to the field line and the 1»5® detector was to look mainly up 
the field line. Since it is generally expected that the majority of 
particles are coining down the field in a precipitation event it is con- 
sistent that the 45° detector always had a higher count rate than the 
135° detector. The 90® detector which crossed over from up looking to 
down looking therefore showed a great variation over a single roll. Of 
the three data sets, that of the 135° detector had the best defined 
maximum which appeared as a narrow peak once per roll. From the preced- 
ing considerations we have made the assumption that this maximum on the 
average corresponds to the lowest pitch angle viewed by the detector in 
its periodic sweep. Using all of the constants previously determined we 
calculated the angle between the magnetic field and the detector look 
direction. (This is the complement of the particle pitch angles) as a 
fimction of ijj for all times of the flight of the rocket during free 
fall. For the correct choice of ij; the peaks of the data set from the 
135° detector should lie on the peaks of the look direction angle and 
the valleys on the valleys. Because of the vagaries of counting statis- 
tics and the time resolution of .00853 second the peaks of the data vary 
from perfect symmetry and exact roll position. To overcome this we 
calculated the cross-correlation of the angle with the data at 0 sec 
lag for various values of . This function is sho^m in Figiire 3.7. 

We chose the value of ^ at the peak of this curve as the best fit 
value and the half width at half maximum as the uncertainty. 

In table 3.3 we list the final values for all the determined 


constants . 
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Table 3.3 

= 1 . 1 ° + . 5 ® 

^ - 1.65° + .10° 

4 )^ = 231.13 + .03° 

V = 78.21° + 2° 

a = 6.049° + . 5 ° 

6 = U8.9° + 1.0° 

w = 1943 + l6°/stc 
n = 8.287 + .035°/sec 


AS X10 


Y degrees 


Figure 3.7, Cross-correlation of the compliiaent of the pitch angle and the count rate for the 
135 degree detector. The angle i|t is shown in figure 2.6. An arbitrary constant was subtracted 
from the actual value of S(iJj) to enhance the variation for plotting. 


CHAPTER IV 


DATA REDUCTION AND EXPERIMENTAL RESULT 
IV. A. Data Overviev 

The data in the nine channels of each of the three spectrometers vas 
recorded simultaneoiosly over a .00853 second collection period. Figure 
U.l is a display of the output of several channels of the three detectors 
over a few seconds of the fli^t. This is simply a display of the un- 
processed count rates. 

Several pieces of information are immediately extractable from these 
figures without any detailed analysis. The strong periodicity that is 
so evident for the 90° and 135° detectors is at the roll rate of the 
rocket, 5*^ rev/second. This indicates a strong spatial dependence in 
the flux of precipitating particles. There is also a distinct mirror 
symmetry about the peaks. Since the detectors s\reep hy the same pitch 
angle twice per roll (going first from high pitch angle to low pitch 
and then from low pitch angle to high pitch angle) this is an indication 
that the spatial structure is ‘basically a pitch angle structure. In this 
particxilar sample the roll rate is not very evident for the 1*5° detector. 
Since it too sweeps out a range of pitch angles this is an indication of 
pitch angle isotropy over the downward hemisphere. There are some 
periods when a weak roll dependence is evident which is the retvim to 
anisotropy. There are also apparent peridicities in the amplitude en- 
velope of all six channels shovni in Figure U.l. We will examine these 
periodicities in greater detail later. 

The munbcr of counts in each chtuincl ns was indicated in Chapter I 
is a measure of the energy spectra. No channel at any time experienced 
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Figure 41 . Unprocessed cour^ rates for several channels A^'D detectors, 
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overflow. In fact coimt rates only exceeded one-half of the maximum in 
channel one dtiring the most intense periods of the fli^t. The spectra 
fell off quite steeply from channel 1. The count rates in channels 
6, T» 8 Slid 9 rarely exceeded 10 counts per accumulation period although 
they had considerably more counts at times than the sample spectra for 
the softer cases in table 1.1 would have lead \is to expect. As will be 
explained, we have excluded the data from channels 8 and 9 from our 
final result because there are good indications that they are contaninate'd 
with proton counts. 

From preliminaiy examination of the data we obtained some idea of 
the direction that our initial ar.alysis should take. There is strong 
time dependent pitch angle structure as well as time dependent energy 
structure. The resolution of spatial dependence from time dependence 
can not be accomplished unambiguously on a moving rocket. The best one 
can hope to do is set some limit on the size of the spatial features and 
ask if they are likely. 

The basic data collection period of .00853 seconds corresponds to a 
spatial distance determined from the rocket horizontal velocity. At the 
top of its trajectory the rocket had a horizontal velocity of .U8 Itra/sec. 
The basic distance then is determined to be about b.l meters. One might 
argue that the Larraor radius is a minimal coherence distance for particles 
in a magnetic field, li.l meters is the Larraor radius of an electron with 
an energy of about 1 Kev. For energies less than this, spatial featiires 
on the order of a Larmor radius are unresolved while for energies greater 
than this such features appear as time variations in the flux. As we 
will see the apparent time structure of the low energy electrons (about 
50 JCev) oJid the high energy electrons (above 600 Kev) is very closely 


85 


i 

? 

I 

i 

i I parallel. The Larmor radiias of a 600 Kev electron is some 4,26 times 
the Larmor radius of a 50 Kev electron and mi^t he expected to give a 
larger coherence length for the 600 Kev electron. Ho such difference is 
apparent in the data. Though hardly a conclusive argianent, this reason- 
ing does take us a hit along the way to accepting the apparent time 
variation in the flux as not arising fhom spatial structure. 
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IV.B. The REP And The FCA 

The BEP of 31 May 1972 began at 6:30 UT as observed by the liometer 
in Andenes. The riometer record is shown in Figure k,2 (Matthews and 
Simons, 1973). 1.0 db of this absojTDtion was due to a solar proton 

event which had begun on 28 May 1972. 

To provide additional criteria for the launching of the rocket we 
were manning the partial reflection facility of the Auroral Observatory 
at Tromso. This experiment provided real time indication of the hei^t 
of ionization in the 0-region. We required that there be significant 
ionization below 75 to before we wo\ild accept an absorption event as an 
BEP. At 6:30 UT there was a marked hardening in the precipitation as 
evidenced by a sudden increase in the absorption of the reflected signal 
from 65 kilometers altitude. This absorption persisted for some five 
minutes, disappeared and then reappeared six minutes later at approxi- 
mately 6:l4 UTi. The ab"orption below 65 km then continued until 6:30 DT. 
The rocket was launched at 06:19 UT and reached apogee at 06:22 UT. From 
the riometer data (Figure it, 2) it can be seen that the absorption had 
already begun to decrease during this period. The returns from the 
partial reflection experiment indicated however, that the detectors were 
taking data during the time of the most sustained hard precipitation. 

As explained in the introduction it seemed quite unlikely that 
protons were the prijnary cause of absorption during BEP’s but it was 
considered essential to include a proton e^qperiment among the complement 
of experiments to be carried aboard the rocket F-34. The proton experi- 
ment was designed at Bell Telephone Laboratories and built , at the 
University of Maryland. ^ the use of coincidence logic between two 
solid state detectors and a magnetic field across the entrance aperture 



the proton detector excluded all electrons from being coxmted (private 
communication L. Lanaerotti, 1972). This allowed us to design our 
electron detector without providing for exclusion of protons. 

Before it was possible to deteimne the actual electron fluxes, it 
was necessary to determine the extent of possible' proton contamination 
of the data. The PCA proton fluxes are shown in Figure 1|.3. These 
measurements are taken from several sources. Three energy points are 
from the Bell Telephone Laboratory proton experiment on the explorer hi 
(E3qp-4l) satelite (private communication L. Lanzerotti, 1973) which was 
located at about 18 r„ in the noon sector. Also shown in Figure 4,3 are 
two points taken from the BTL-proton experiment on P-34. We show 
in addition the average flux calculated from channel 8 of our electron 
detector assuming the counts were all due to protrons. 

The equivalent energy thresholds for protons detected in the electron 
spectrometer were quite different from the thresholds for electrons (see 
section I-F) , We show the equivalent thresholds for channels one through 
nine in table 1.6. The proton fluxes are insignificant when compared to 
the electron fluxes for channels one through six. The average flux of 
electrons in channel 7 determined for the entire flight is possibly 
eont«uninated by as much as 15%. We have determined this value hy extra- 
polating between the points from the proton experiments shown in Figure 
4,3. It is obvious from Figure 4.3 that all of the counts in channel 8 
and 9 may be accorinted for hy the protons, though an extrapolation be- 
tween the Exp-4l points indicates that about 7^^ of the counts in channel 
8 are attributable to electrons. In table 4.1 are shown the predicted 
average count rates in channel 7j 8 and 9 based upon Explorer 4l satellite 
data (L. Lanzerotti, private communication) Compared to these are the 
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Figure 4.3, Proton fluxes due to Polar Cap Absorption as determined 
from several sources. BTL-Bell Telephone Laboratory proton ^periraent 
on F-34, EXP-41“Bell Telephone Laboratory proton detector on satelite 
Explorer 41* UOM-Unlversity of Maryland electron experiment of F-34. 
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IV. C, RemovinK The Roll Dependence. 

The time structure so obvious in the raw data indicated that un- 
folding the roll dependence was a reasonable first step. As the detector . 
sampled the complete pitch angle distribution twice for each roll of the 
rocket we decided to take one roll as the basic data period assuring two 
passes at each pitch angle. If the time variation in the flux is slower 
than tne roll rate then our averages are well representative of the 
actual time structure in the pitch an^e. If the time variations in the 
flux are greater than the roll rate then the results are a frequency- 
filtered component of the true time variation. ¥e will examine this point 
more fully in section U.F. 

As the detector swept through -various pitch angles during a data 
collection period the central look direction vector of the collimator, 

Vjjj , pointed to different pitch angle as a function of time. The curve 
of small dots in Figure U.h shows the pitch angles covered by the vector 

for each of the three detectors. In this curve three dots delineate 

each data collection period; one for the initial pitch angle, one for 

the pitch angle half -way through the collection period and one for the 

final pitch angle. A careful examination of the figure shows that the 

points are not equally spaced in pitch angle. This unequal spacing in- 

% 

dicates that the width of the pitch angle field of -view is different for 
different data collection periods and this width is a function of 
position along these curves. It can also be determined from Figure 4.4 
that although the detectors never looked at the same pitch angle simul- 
taneously they did none-the-less have overlapping pitch_ angle coverage. 

During a single sampling period particles enter the collimator from 
all pitch angles within the field of vici;. Each pitch angle contributes 
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flux in proportion to the area exposed to it and the length of time the 
area is exposed. If we let F he the total flux of particles that 
enters the detector during a period At then 


t+At 

P = I I [G(a»t)J(ast)3 da at 

t a 


(4.1) 


where G(a»t) is a geometrical term determining the pitch angle weight- 
ing or response of the detector at time t . J(a,t) is the integral 
flux at pitch angle a and time t , In terms of the geometrical 
quantities defined in Figure 4.5 we may also write 


t+At 


-niH 




J ( a, t ) r^rgdr^argd<{)^a(j>2at 


(4.8) 


2 ^1 -"2 


with = r^ + y^ - 2r^rg cos(4>j^-4>2) ^ cosY « h/D . Ihe vector 

A. A 

passes through the differential areas and Og malting an angle 

a with B 9 the magnetic field. The angle Y is the angle Between the 
magnetic field and the look direction vector* * of the collimator. 

jSiS should he obvious from equation (4.2) and the considerations of 
Appendix A-4 the pitch angle response discussed here is not the same 
function as the angular response of the collimator discussed in section 
I-D and is in fact a more complex function. 

By taking the dot product of ip and B we can obtain a relation 
between the pitch angle a and the variables of integration in equation 
4.2. 
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(r^siti4'j^~rgSin4>g) sinY + h cosY *= D cosa 


This equation may be used to replace one of the four variables .r^ , 

Tg 9 4'i i or ({>2 in equation ^>.2 with a so that a becomes a variable 
of integration as in equation U*l, Mahing this substitution we obtain 
the result 


t*lAt g ft 

^'1 I I i \ ^ 

” 


t a Yg<l»i4»2 


where we may immediately identify the function G(a,t) as 


G(a.t) = | I I 


rg 4.g 


With the time dependence entering through the function 


= r^(a,Y(t)»rg,(j)^,4i2) 


(4.7) 


For any given data aquisition period we may deteimine GCajt) by 
numerically integrating equation 4.6 for a fixed value of a and Y * 

We have instead determined the fiaiction G(ajY) » (the time averaged 
value of G(a,t) over the counting period), which we will refer to as the 
pitch angle response function, G(a,Y) is given by 


t+At 


[a,Y) - ^ I G(a,t)dt 


(4.8) 
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vhich we associated with the value Y given "by 


t+dt 

At j Y(+)at (U.9) 

t 

in Figure U.6 we display G{a,Y) for several values of Y . As we 
used G(a,y) as a wei^ting function over each counting period it was 
always normalized so that it had a maximum vadue of one for each period. 

With the pitch angle response function, G(«,y) * the pitch angle 
dependence was \mfolded from the data for each of the three detectors. 
Since this involves talcing data from approximately 22 consecutive data 
points the result at hest represents a time average over O.I85 seconds 
and a pitch angle smoothing over 15 degree intervals. 

Khe measured flux in any collection interval, i , which is given 
in terms of the number of counts, , the geometrical factor, G , the 
collection time At and the efficiency E hy the relation 



(4.10) 


is assumed to first order to come from all directions in the acceptance 
view of the collimator. As a statistical weight for the flux at 
pitch angle a , we took the pitch angle response at a , 6(a,Y)« ®ie 

average flux at a given pitch angle during one roll is given hy the 
expression 



(ll.ll) 


22 _ 

I F G(ci*Y ) 

-/ ^ i=l ^ ^ 

J(a) = 

I G(a,7.. ) 

• *1 «L 

1=1 


It follows that the iincertainty is given by 


^2 I mT T 

« yi=l ^ \(GxAtxE)^/ 

®j(a) “ 22 _ 

I G(a,Y ) 

i=l 


(4.12) 


We have asstmied 0„. = v4l7 and that there is normal propagation of error. 

1 

For the tvrenty-two points in a single roll of the rocket, G(a,y^) is 
Iqrpically zero for a given cc at all but 6 values of • There are 
therefore only about six terms in each of the above siims for a specific 
pitch angle a , 

For each channel of each detector we ccuid then produce pitch angle 
distributions of the type demonstrated by Figure h.J. There are two data 
sets shovm in each diagram in Figure 4.7- They are for corresponding 
energy channels for the 45 degree and 90 degree detectors. For much of 
the precession period of the rocket these detectors had overlapping pitch 
angle coverage during any one roll of the rocket. The detectors did not 
look at the same pitch angles simultaneously nor did the corresponding 
energy channels have the same response functions. We woiild not therefore 
e35>ect them to record the same integrated flux at the same pitch angle. 
This accounts for the apparent inconsistency in the overlapping regions 
in Figure 4,7. We do ei^ject that the representative energy spectra ob- 
tained at the overlapping pitch angles be equal within the limits of 
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Figure 4 . 7 . Typical pitch Af^LE DisTRiEirncHs DETERi'iifJED oy'er one roll of "He roqcet. 
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uncertainty. Pitch angle distributions were produced for every roll of 
the Socket from time 80 seconds to time 300 seconds into the flight. We 
will discuss these distributions in section 4.E. 

The pitch angle distributions enabled us to obtain an overall pic- 
ttire of the intensity variations on the time scale of one roll of the 
rocket. ¥e defined the average unidirectional flux, , over the pitch 
an^es covered by each detector. J is defined by the relation 


j * 


2ir am. 

J I J(a)sina dad<}» 
0 omin 


Sir oma 

J J sina dad4> 
0 omin 


(4.13) 


omin and amax aj.‘e the pitch angle limits of the pitch angle 

~2 -1 -1 

coverage for a given detector. has units of cm x sec x ster 

Figure 4.8a shows for the first half of the fli^t for two channels 

of the 45° and 90 ° detectors. At 80 seconds the rocket had reached an 
altitude of 90 km so that even in the lowest energy channel there was 
little if any atmospheric energy degradation of the particles reaching 
the detector, A 4o kev electron will penetrate to about 8T km leaving 
most of its energy in the last 4 or 5 kilometers (Maeda, 1965). The 
fairly constant flux level for the first 15 seconds of data shorm in 
Figure 4.8a is a feature of the precipitating particles therefore and not 
on atmospheric effect. At about 107 seconds there is a marked increase 
in the level of particle fliix characterized by a steep increase in the 
counting rate. This steep increase in the counting rate occurred simu- 
ltaneously in channels 1 through 7 of all three detectors to within the 
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Figure 4.8b* The average unidirectional flux, J • in channels 6 through 9 o£ the 45* 
detector (det-2) . '* 







102 



time resolution of .185 sec {one roll), The fractional increase in flux 
is much greater in the higher energy channels and can he seen in Figure 
U.8h which shows the average unidirectional flux for channels 6 through 
9 of the 1(5® detector. As is somewhat more evident from the higher 
energy channels there is a decrease in flux at 130 seconds and then an- 
other increase continuing until a rather abrupt drop-off at 170 seconds, 
■We refer to these periods as htursts 1 and 2 . There is another such hurst 
beginning at about 230 seconds (Fig. h.8c) continxiing to 265 seconds and 
another from 278 seconds to 288 seconds. We call these bursts 3 and 4 . 
The bursts are characterized by pulsations in intensity over which the 
flux changes by a factor of two in the lowest energy channel and by as 
much as ten in channels 6 and 7* Superimposed on the flux for the entire 
fli^t are short period pulsations which involve both large and small 
changes in the flux. We will refer to these as the fine structure . We 
will return to an analysis of these various features of the time struc- 
ture in section IV-E and IV-F, 
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IV.D. 


R,(B) is the response function of the ith channels, as we have discussed 

*L 

in Chapter U, Q is the geometrical factor, and At is the accumulation 

time. The ith accumulation will yield a measured count N. . To fit 

im 

an analytical foim to F(E) we will define the chi-square parameter 


« T (N. -M.)‘ 

= I 

1=1 im 


(h,13) 


which will permit us to determine a least square fit to the measured 

data. The uncertainty for N. will he assumed to he from 

im im 

p 

Poisson statistics. t|i must he minimized with respect tp the constant 
parameter that define the function F(E) . 


The hroad response functions of the energy analyzers (see Figure 
2.16) were capable of distinguishing only large features in the energy 
spectra. It is the common practice when looking for hroad features to 
attempt to fit exponential spectra (as we discussed in Chapter II), power 
law or some similar general form that can give parameters that character- 
ize the energy hardness and flux level. Assuming a spectral form, F(E) , 
which over any accumulation period is the average unidirectional differ- 
ential flux over the aperature, the number of counts in ‘^■he ith channel 
of the analyzer is given by 


Emax 

= G-At*x| R^(E) F(E) dE 
0 






Both poorer law and exponential law fits to the spectra were attempted. 
They proved generally Inadequate in representing the spectrum over the 
full range of energy. The exponential law was generally "better than the 
power law. Figure U.9a, a best fit exponential law is shown by the 
solid line. The corresponding point snectman is also shown with the 
appropriate error bars. We explain a little later how the point spectrum 
is determined. The exponential law consistently gave too large a flux 
in the lowest energy channels and too small a flux in ■Mie_ higher energy 
channels. The power law was generally inadequate across the entire 
spectrum. A very good fit to the data could always be obtained by 
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combining the exponential and power laws into a four parameter fit of 
the form 


F(E) » E”* + Jg Qxp {- E/Eq} (ij.l9) 

but this complicated the iterative fitting procedures lased on the computer 
to such an extent as to require a prohibitive amotmt of computational 
time to process the data from the entire flight. In an attempt to pro- 
duce a better representation in reasonable computation time we developed 
an ad-hoc procedure for defining the spectrum which reproduced essentially 
the same spectrum as the four parameter fit given by equation (lj.l9). 
tChis is the point spectrum shown in Figures ij.9a and b and refered to 
above. This procedure guaranteed that the spectrum gave the correct 
count rate in each channel so that the error for the assigned flux was 
taken from the counting statistics. The error in the energy was taken to 
be the full width at half maximum of the response function of the channel 
folded with the resulting spectrum. 

We assume that the spectrum in two adjacent channels may be re- 
presented by a power law. 


T = E~* (U.20) 

so that the counts in the two consecutive channels are given by 

Eoiax 

Nj « GAt I R^(E) dE - - (lt.21) 

0 


and 
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Emax 

0 


(U.22) 


The ratio of coiaits 



Emax 

R^(E) E”'^ dE 

0 

Emax 

I R..,(E) E*** dE 
0 ^ 


( 4 . 23 ) 


is strictly a function of the hardness and thereby automatically deter- 
mines X . Computationally a table of ratios may be determined as a 
function of x and the values of x sinqply read from the table once the 
ratio is knomi. is determined by 



( 4 . 24 ) 


once X is Imcn'm, With the values of and x determined, the 
average energy of a count in channel i is then given by the relation 




E E”* dE 


E”^ dE 


( 2 . 25 ) 


The value of the spectrum is evaluated at E. and the spectral 


point then chosen as 




Spectra deterioined in this fashion are displayed in Figures h.9a and h. 
Ihe power law determined over each channel is drawn as a dashed line 
throu^ the points. Figure 1^.9h is the average spectra determined from 
the data over the entire flight. 
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IV. E. ignergy And Pitch Anfcle Structttre. 

The reduced data from the REP experiment contains a tremendous 
amount of detail which woxild prohahly justify further months or years 
of labor to more fully elaborate and understand. Several broad features 
have emerged which we feel may fairly well illuminate the physics of the 
REP phenomenon. 

Figure ^t.lO shows a time display of the average integral unidir- 
rectional flux of precipitating particles shown in parallel with the 
average energy of these particles. In order to improve upon the time 
resolution of the data in Figures U.8 these averages have been determined 
over a period of sli^tly more than one-half roll of the rocket .095 
seconds. The time structure which emerges is much the same as that in 
Figures U.8. There appears to be two types of precipitation; that during 
the bursts 1, 2, 3 and 4, and that during the interburst periods. The 
bursts are made up of short-bursts which have a fine structure super- 
imposed upon them. The interburst periods have the fine structure 
pulsations which are superimposed on small bursts similar to the sub- 
structure of the large bursts. The burst and pulsing behavior in the 
integral flux is paralleled by a consistent behavior in the pitch angle 
distributions. In Figure 4.11a we show the pitch angle distribution at 
the peak of a large burst for several energies. In Figure 4,11b we show 
the pitch angle distribution at the valley between two pulsation peaks. 
TOiere is pitch angle isotropy at the peaks and a loss cone anisotropy at 
the valleys. The loss cone is more pronounced for the higher energy 
particles. The consistency of this behavior is demonstrated by the time 
plots of pitch angle for the precipitating particles shown in Figures 
4.12a, b, c, and d. The Figures 4.12 are for channels 3 and 4 of the 43- 
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Figure 4 .lie. The pitch angle distribution (interpolated to three energies and normalized 
to approximately the same peak value) shown at the peak of a pulsation. 
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Pi 5 sure 4.11b. The pitch angle distribution (interpolated to three energies and normalized to 
approximately the same peak value) shown at the valley between pulsations. 












degree detei^tor taken over the periods of burst 1 and 3. The pitch 
angle distribution at the top of the bursts is more nearly isotropic the 
higher the flux level, but even at the lower levels there is a definite 
tendency toward isotropy as the flux increases. This behavior seems to 
be qualitatively the same for the large burstsj the small bursts, and even 
the fine structure pulsations. 

Figure it. 10 shows that the average energy increases and decreases 
ri^t along with the flux. This is of course a consequence of the more 
pronounced increases in the flux of the higher energy particles which 
may be seen in Figure 4.8. To within the time resolution of Figure 4.10, 
,095 seconds, the energy and the integral flux have a cross correlation 
which peaks strongly at zero. The change in the energy spectrum which 
accompanies the time variation in average energy is portrayed in Figures 
4.13a and b where the spectra have been determined at adjacent peaks and 
valleys. The log-log scale in that figure may be misleading but a close 
examination shcn^s that the proportional enhancement in the fliix level is 
an increasing function of the ener^, 

A very clear picture of the precipitating electrons emerges from 
Figiures 4.10 through 4.13. Particles appear in the loss cone simultan- 
eously with an overall increased fliuc and an increased hardness. The 
increased hardness is brought about mainly by the enhancement of the 
high energy tail of the distiibution. These particles out in the high 
energy tail are of coiirse responsible for the event having been a REP 
in the first place, and we will pursue their origin in Chapter 5* We 
point out for exphasis at this point that most of the clear and positive 
results obtained fjpom om* data analysis are outlined in Figures 4.10 
through 4,13. 
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IV. F» ^Pime Series Analysis 

The level of fluctuations which is so evident in Figures U .8 and 
U.IO suggests the possibility of strong freq.uency dependence in the REP. 
There are various forms of anagneto spheric activity that have distinctive 
frequency characteristics which mi^t be expected to appear in associated 
particle precipitation. It therefore seemed incumbent upon us to check 
the data systematically for frequency information. 

As indicated previously the basic data accumulation period was 
,00853 seconds. This time period sets a- theoretical limit on the fre- 
quency content of the data of about 58 hertz, (see below), While this 
is a very low frequency conpared to the possible electron plasma and 
cyclotron frequencies that one expects, it certainly includes the range 
of geomagnetic aaicropulsations and ULP cyclotron waves that have been 
observed in the magnetosphere . Electron bounce times for Kev electrons 
are on the order of a second as well. 

The rocket roll frequency is the dominant frequency present in the 
count rates from each of the detectors. This frequency is present be- 
cause of the pitch angle anisotropy that existed to at least some small 
degree at all times during the fli^t. This periodicity is special. 

Any other special periodicity that exists in the pitch angle distribution 
will also be present in the raw count rates i|ust as the roll periodicity 
is present. This sort of frequency information is indistinguishable 
from the overall fluctuations in the flux that may he present at all 
pitch angles. This is somewhat unsatisfying in that we would like to 
resolve the question as to whether the variations in flux jLevels that are 
present in Figures U .6 are trxaly the basic time structure in the event 
or a strongly filtered component. The averaging time period for the data 
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in Figure U.6 is .095 seconds which yields an upper freq,uency cufcoff of 
about 5Hz (see below). We will demonstrate by a simple example how the 
hi^er frequencies are filtered by the averaging. Shis filtering does 
have the advantage of decreasing the effects of aliasing. 

We suppose that the true time behavior of the flux may be represented 
by a function f(t) with a Fourier transform f(o) , which is related 
to the spectral density S(oa) by the relation 

S(to) = ||f(w)|^ (4.26) 

0? is a period long with respect to all inverse frequencies of 
interest. When we sample the function f(t) by time averaging over a 
period 1? we actually create a new function f (t) defined by 

t+T/2 

fg^(t) « i I f(tO dt‘ (4.27) 

t-T/2 

which is sa3i5>led at evenly spaced points t. . We trivially determine 

the Fourier transform of f (t) by substituting the Fourier transform 
. ^ 

of f(t) into equation (4.27) interchanging the order of inter- 
gration. 


f Jt) 
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Carrying out the integration we find that 


(U.28) 





f (t) = 


sin(<ifl^/2) .“ioJfc ,.. 

“T<*72) ® “ 


Frcffli which we identify 


■p (,^\ - sin(t»fl'/2) 


and the power spectral density 


S „(tii) = S(u) 


sin^(tJI?/2) 

(ftST/2)^ 


Those freq.uencies which satisfy the relation 


= rnr 


are totally absent in S (w) . The half power point is at the frequency 

av 

*\/2 relation 


^2 ^^1/2^^ ^ 


with the power dropping off veiy rapidly h^ond this frequency. (See 
Figure U.lU) 

In the digital technique applied to determine the power spectral 
density we use a finite set of equally spaced data points. The useful 
frequency region will he from 0 to ir/T radians per "second since 
ftequenoies higher than tt/T will he folded into the lower frequency 



range. The upper cut-off frequency is called the Nyq.uist frequency 


las 


w = ir/T (U.3U) 
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The frequencies w , higher then the Hyquist frequency which are 

8f * 

folded into frequencies, to , helow the Hyquist frequency, are called 
aliased frequencies and they satisfy the relation 


to 2nto + to 
a c — 


(U.35) 


with n being an integer (Eendat and Piersol, 1966). In our case alias- 
ing is much less significant in that to is si j ghtly greater than the 

o 

half power frequency, bhat all frequencies above to^ atre 

fairly well filtered. 

The frequency analysis was performed using standard digital tech- 
niques which may be foimd in Bendat and Piersol (1966), The definitions 
are repeated here for coisqjleteness . 


T is the basic data sampling interval and the deviation of the 

nth point in the series x^^ is given by 
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X - X 
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(U.36) 


with X being the average value of x^ and x^ being the value of 

x(t) at t =! nT . R is defined to be the autocorrelation coefficient 

m 

at t = mY and is given by the relation 


H-m ", n n+1 
n-j. 


We may then determine the raw estimate of the power spectral 
densilqr S^(f) . 
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Ih which we take f as the frequency in hertz so that f = w /2w 

c c 

M is the maximum value of m where m is referred to as the number of 


lags, f in equation 4.34 is evaluated at even multiples of 


, ^ 


our estimates the raw estimates have "been smoothed with ’Harming' 
formulas 


S(f ) = .23 S (f ,) ^ .54 S (f ) + .23 S (f 
n n—-L r n r 


S(f^) - .54 S^(f^) + .45 S^(f^) 


(4.39 


p(f„) = .5!. S^( V + M 


The standard error of S(f ) is given hy the relation 

n 




(4.40 


so that the maximum number of lags, M , must he chosen small with re- 


spect to the number of pointy H 
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She spectral analysis vas performed on hoth the averaged data and 
the ravr data. She averaged data had a time interval of .0938 seconds be- 
tween points and therefore a Nyquist frequency of 5.33 Hz. Shis is a bit 
less than the rocket roll frequency of 5* ^^3 Hz. The Nyquist frequency 
for the raw data set is 58.8 Hz. Figures U.15 a * b show an exan^le of 
the power spectral estimates over the first third of the fli^t for two 
energy channels from the averaged data set. There are four peaks which 
show quite distinctly in these figures and are also present in the 
spectral estimates for all of the channels one through seven. These are 
the peaks at .2, *SU, .6U and .906 Hz. In addition there are peaks at 
.^53, .800 and I.U9 Hz which are more distinct the higher the energy of 
the channel. The standard error is indicated on the figure. As might 
be eaqoected* the roll frequency has aliased into the spectral estimates 
at 5.2 Hz. As a result of the frequency filtering though its contribu- 
tion to the power spectral estimate has been reduced by 62 per cent. All 
of the piominent peaks are sufficiently below the half power frequency 
so as to not warrant any co37rection. 

The most significant feature of the power spectra in Figures 
k.ll*a and U.lkb is the lack of any really dominant frequency. The fine 
structure time variation certainly appears to have a periodic structure. 
The power spectra indicates that there are in fact several low fre- 
quencies sti’ongly present in this quasic-periodic behavior. 

Figures 4. 15a and b show examples of the power spectral estimates 
over the first third of the flight for two energy channels from the raw 
data. The largest peaks are from the roll variation which shows up out 
to as much as 5 harmonics. These very strong features tend to swamp any 
features too close to thesi. Just resolved from the first harmonic of the 
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Figure 4* 15a. Power spectral estimates ovef averaged data set for channel 3 
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roll is a pealc at 6.17 Ha which also appears at the second harmonic 
frequency of 12.34 Hz. 6.17 Hz is present in channels one through seven 
from the beginning to the end of the flight. There are evidently many 
peaks in the spectra which are statistically significant, thotigh none of 
them dominate (with the exception of the roll frequency). These results 
leave us with more of an indication of what is not present in the data 
than any positive indication that we had hoped might emerge. 

Using teehniq.ues similar to the power spectral estimates, cross- 
correlations were determined between the several energy channels. The 
expectation is that a precipitation mechanism acting over a finite region 
of space, sufficiently removed from the rocket, will lead to arrival 
time dispersion in the particle counts at different energies. As with 
the above results the information was of a negative rather than a posi- 
tive nature. To within a time resolution of .00853 seconds the particle 
fluxes in ail channels were most strongly correlated at a log time of 
0 seconds. If 50 Kev and 600 Kev particles were simultaneously injected 
into the loss cone at the equator on an L shell of 6 , they would arrive 
at the foot of the field line with a time separation of about ,196 
seconds. This is a time which is resolvable in our data. We have the 
problem though that the roll correlation my be so strong and broad as 
to swamp all other correlations that lay near hy. The energy width of 
the various channels is also broad and would tend to broaden and decrease 
any true correlation of this sort. Lastly the finite width of any inter- 
action region about the equator would also lead to a smearing of such a 
dispersion effect. We are left to conclude that the absence of such 
dispersion related correlations is not in itself definitive. Our time 
resolution is such though that if we were to take these results as 
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Figure 4.1Cb. Power spectral estimctes raw data set channel 4 for lOjQDO data pointsi 






CHAPTER V 


SOME THEORETICAL CONSIDERATIONS ON REP’S 
V.A. The REP And The Trapped Particle Population 

As has been indicated in the introduction, the REP is quite uncommon 
when compared to magnetic substorms (Bailey, 1968). If REP’s are truly 
a substorm related phenomenon then some other features of the magneto- 
sphere must determine when substorms will lead to REP’s. Numerous papers 
have been published detailing the enhancements of trapped particles in 
the magnetosphere during and following magnetic substorms (see for instance 
Pfitzer and Winckler, 1968). The reports based upon data taken by 
synchronous satellites which have been located on the L = 6 field lines 
(Mozer, 1972; Bogott and Hozer, 1974) are of particular interest in this 
respect. Rosenberg et al (1972) compared the precipitated flux with the 
trapped particle flux at synchronous orbit for the REP events of 25 
August 1967 and 11 August 1967. They determined that there were marked 
enhancements of the trapped particle population prior to the REP's 
probably due to substorm associated injection of particles on the night 
side of the magnetosphere. While the 25 August event xtas accompanied 
by significant enhancement of trapped relativistic electrons, the 11 
August event was not. They concluded that there were at least two dis- 
tinct subclasses of REP’s. The best e-folding energy which would be 
fitted to our average fluxes was about 40 Kev. This corresponds very 
well with the hard component of the precipitated flux of the 11 August 
event (Rosenberg et al 1972). The e-folding energy of the hard component 
of the 25 August event was significantly greater than "this. Though it 
is not conclusive the data does support that the REP of 31 May 1972 was 
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of the same class as the BEP of 11 August I 967 and therefore might 
not have been associated with a significant enhancement in the trapped 
relativistic particle population. In our following considerations we 
have not made any attempt to draw any distinction between the two types 
of EEP's but rather looked at the more general problem of how might 
relativistic electrons be produced in a post-substorm magnetosphere. 

Many authors have discussed electron-cyclotron wave gyroresonant 
interactions (for example Kennel, I966; Kennel and Petschek lS’66j Kennel 
and Engelmann I966; Lerche I968; Coroniti and Kennel 19 T 0 ) and it is 
generally accepted that these interactions are a major cause of diffusion 
of electrons into the magnetospheric loss cone. As described by Kennel 
and Petschek (1966) and many authors since, electron whistler interactions, 
imder certain conditions, establish an upper limit to the stably trapped 
fluxes. The atmosphere at the foot of each field line establishes and 
maintains a loss cone anisotropy. Particles diffusing down into the 
loss cone tend to spread the anisotroxy to pitch angles above the atmos- 
pheric loss cone. Whistler waves are unstable in the presence of such an 
anisotropy and will grow forcing particles down into the loss cone tending 
to remove the anisotropy. A complicated nonlinear equilibrium is 
established in which wave loss and growth and particle loss and diffusion 
all balance each other. A completely self consistent and general mathe- 
matical solution to such a problem has not been carried out. Kennel and 
Petschek (1966) attempted a solution but were forced to make many assump- 
tions about the nature of sources and the form of the pitch angle 
anisotropy and therefore intended their model to he only semiquantitative. 
A long term study by Mozer (1972) of the trapped particle population on 
h - 6 showed that electron fluxes rarely exceeded the limit predicted 
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"by Kennel and Petschek (I966) and never exceeded that limit by more than 
a factor of two. The integral flux of electrons (E>50 Kev) that was ob- 
served during the REP of 31 Ifey 1972 was at times at or slightly above 
the Kennel and Petschak stable limit. The particles observed at I50 km 
are mainly in a region of pitch angle space which corresponds to the 
equatorial loss cone and therefore normally exhibit lower fluxes than 
the trapped flux. It seemed reasonable following these considerations to 
examine the cyclotron resonant diffusion theory to determine if the in- 
clusion of relativistic terms mi^t not point the way to a possible pro- 
duction mechanism for EEPs. A number of very interesting and possibly 
important points have emerged from oirr investigation, 

V.B. The ^-Jbistler Dispersion Relations 

The standard treatment (as exemplified hy the above references) 
of the parallel propagating right hand cyclotron wave mode (referred to 
also as R-H mode and whistler mode) assmes a cold backgrovmd plasma, 
the dielectric porperties of which deteannine the wave dispersion relations. 
We say the cold plasma supports the plasma -vraves. The interactions of 
these 'cold plasma waves* with the distribution of highly energetic 
electrons aire the wave particle interactions of interest. It is assumed 
that the effective number density of the ’superthermal ' electrons is so 
low as to not affect the dielectric properties of the 'background' 
plasma. The complete dispersion relation for the RH-mode waves may be 


written 




u 


; i 

I 1 


If 






135 


(Krall and Trivelpiece, 1973) with u) (u , ) and (H . ) being 
^ pe pi ce Cl 

the electron (ion) plasma frequency and the electron (ion) gyrofrequency 
respectively. These frequencies are given by the relations 


rt e _ 

. . 2 _ a , „ _ eB 

u and j2 „ 

pa ca m^c 


5.2 


with a standing for the ion species, the ion mass and n^ the 

number density of species a * (We are using cgs gaussian units and 
will do so throughout this chapter and associated appendices). 

For a neutral hydrogen plasma one has n^ - n^ and = 5.^ x 

..li 

10 so that equation 5.1a may be written 
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5.1b 


For virtually all discussions of electron-whistler mode interactions 
in the magnetosphere the articles in the literature contain an approxi- 
mate form of this dispersion relations. The assuiiq>tion is made that 


ft - < w « 

Cl ce pe 


so that the dispersion relation may be written 
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When condition 5.3 is not met equations 5* lb and 5.1c yield greatly 


different results. In fact at the equator on the L = 6. field line 
there are times when condition 5.3 is not true. Typical quiet time 
densities may run as high as 5/cm'^ (Chappol et al, 1970) while dtiring 


Q 

distiirbed periods following substorms the density may be less than .1/cm 
(Williams and Lyons, 197^)* Using a dipole field approximation the field 
at the equator is given by 


B = gauss 

o L 


5.4 


This yields a value for £2 of 26 . rad/sec. For n = 5/cm , 

ce o 

3 

0 ) St 126 rad/sec and condition 5.3 can be met. For n^ = 0.1/cm , 
pe o 

(I) l8 rad/sec so that condition 5.3 cannot be satisfied. In Figure 
pe 

5.1 we show the exact dispersion relation and the approximate dispersion 

relation for ease of W « l8 rad/sec. For frequencies, (O, below 

pe 

0.1 £2^^ the two dispersion relations are in fairly good agreement. 

The range of frequencies that are dealt with in a theoretical treat- 
ment of the wave-particle interactions depends upon the energies of the 
particles of interest. These frequencies are determined through the 
resonance condition which is obtained in appendix B (see equation £.42). 
The resonance condition may be written 




K 


5-5 


in which v_ is the particle velocity parallel to the magnetic field 

and Y is the relativisitc mass correction, to and K are related 

through the dispersion relation. We will consider particles in the tens 

3 

of Kev range to the near Mev range. With n^ = 0.1/ cm equation 5*5 
combined with 5.1b tells us that a 60 Kev electron is resonant with a 
wave with to > . 5 J2 . This is well out of the range of validity of 

C6 

the approximate dispersion relation. 
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Vrom the above we are lead to conclude that during disturbed times 
when particle densities were low the full EH-mode dispersion relation must 
be used to model cyclotron wave-electron interactions. It is also impor- 
tant to point out that any theoretical predictions made using the approxi- 
mate dispersion relation are at least questionable at times of low 
density exterior to the plasmapause. 

V.G. The Quasilinear Diffusion Equations 

In the following we will use the full cold particle dispersion 
relation (eq. 5«lb) and the fully relativistic guasilinear diffusion 
equations, which we have derived in appendix B, to examine the electron- 
whistler interactions in a low density case as described above. In 
deriving the following equations we have assumed that there is a region 
of sufficient uniformity about the equator ro treat the magnetic field 
and plasma as being equivalent to that in a uniform infinite plasma. 

This is a standard approximation the justification of which is to be 
found in the article by Kennel, I 966 and Kennel and Petschek, I 966 . 

The set of quasilinear equations (see appendix B) are 
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in which. Wj(k) is the linear growth rate of waves of frequency w(fc) , 

E^ is the spectral energy density of these waves, “pgh plasma 

frequency associated with the number density of supertbermal electrons, 

f^^ is the distribution function of superthermal electrons, dJlj^ is 

a small solid angle about the parallel propagation vector of waves 

k and Pj^, pn, v, and Vh are the racraientum and velocity of the particles 

perpendicular and parallel to the magnetic field, f^^ is a function of 

time and momentum but is defined in such a way as to be spatially uniform 

(see appendix B). Azimuthal symmetry is assumed about the magnetic field. 

9f /3t is the quasilinear change in the distribution function due to 
00 

wave particle interactions. Equation B .50 is essentially a diffusion 
equation in momentum space. We have derived this equation by following 
the basic procedure of Drummond and Pines (1965), although the same 
results should be obtainable using a Foklie r-Planck approach. To the 
extent that the approximations are appropriate to this set of equations 
should describe both the energy and pitch angle diffusion of the electrons 
due to interactions with whistler mode waves. It is a trivial exercise 


to show that in the limit of nonrelativistic particles and pure pitch 
ancle diffusion equation B .50 is identical to that given hy Kennel, 1966. 

For our later considerations we will modify equation 5.6 to the 
equivalent differential form 

5.6E 


and introduce a convective terii which will allow wave energy to be 
removed from the equatorial ijiteraction region at the wave group velocity. 
We write this as 


aE- V 

dt L 


5 . 6 c 


in which v is the wave group velocity and L is a length character- 
S 

izing the size of the interaction region. 

As formulated above f^^ is a function of and p« where Pn 

extends from -«> to The dispersion relation 5* lb has two branches 

one for positive going waves and one for negative going waves, (Krall and 
Trivelpiece, 1973). We will make some symmetry arguments which will 
limit the domain in consideration as well as simplify some of our 
further discussion. We argue that there is no preferred direccion para- 
llel to the magnetic field so that as many particles move one way as the 
other. This is also true for waves. This leads to the following con- 


clusions: 
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If “we can reduce our considerations properly to the development 

of f for positive Pn we will automatically know it for - pn. 
eo 

This should he true for the vave spectrum as well. Knowing the spectral 
energy density of positive going waves is equivalent to knowing the 
wave spectrum for negative going waves. 

Consider equation P.U 5 above for the linear growth rate. This 
expression must he evaluated with the resonance condition, equation 5.5» 
at each value of pj^^. Eqiiation 5*5 niay also he written in terms of the 
resor^iint momentum 
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As YU} may he larger than SI (provided that y is svifficiently large) 
there will he both positive and negative value p^. In the non— 
relativistic limit the resonant momentum is always negative since O) 
is always less than Q (for cyclotron waves). We will restrict the 
range of p„ to only positive values in equation B.^5 "by breaking 
equation 5*9 into tvro conditions 
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vlth Pj, alvays greater than zero. Equation B.liS then is reuritten as 
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in which we have used the symmetry conditions 5.8. We simply consider 
what happens to the positive going particles and assume the same must 
happen to the negative going ones. 

V.D. The Implications of Pitch-Angle Anisotropy 

At this juncture we will attempt to motivate the rest of our dis- 
cussion hy considering a particular class of distribution functions, 

f (p,» Pii) and the implications of equation 5.H for the linear growth 
eo “ 

rate with respect to this class of functions. It has been noted experi- 
mentally that substorm associated injections of energetic particles 
tend to be at flat pitch angles (Bogott and Mozer, 197l»; Pfitzer and 
Winckler, 1968). The atmosphere at the foot of the field lines main- 
tains a loss cone anisotropy even as other magnetospheric interactions 
work toward removing it. We will consider then a f^^ which is a 
monotonically decreasing function of energy but with a flat pitch angle 
distribution. Such a distribution may be written as 


"l f ^ 

in which Apj_ and Ap„ are parameters which determine the relative 
partition of energy perpendicular and parallel to the magnetic field 


and is an effective particle number density. It is noted that 
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f is not strictly speaking a true loss cone distribution in that f 

0O 00 

is not zero for pj^ equal to zero. This may be corrected by multiply- 
ing fg^ by (sina) where a is the pitch angle and £ is an 
integer. This only complicates the following analysis and does not 
change the conclusion we will arrive at. We will return to the loss 
cone problem later in our considerations. 

f has the property that it can be made anisotropic perpendicular 
to the field line (a flat distribution), isotropic and anisotropic 
parallel to the field line (an oblong distribution) . This is effected 
by controlling the ratio of Apj^ to Ap,|. We will obtain flat distri- 
butions by always having Ap, > Ap„. f^^ also fulfills the necessary 
condition that it satisfy the zeroth order linearized Vlasov equation 
which is equation B.5 in Appendix B with the right hand side set equal 
to zero. 

Expression 5«ll for the linear growth rate represents a balance 
of energy. The waves interact with particles gaining energy frcen some 
and growing and losing energy to others and damping. The contribution 
at each value of pj_ is determined by the expression under the integral 
sign. If the expression is negative then the waves are damping and 
therefore losing energy, if the expression is positive. t^e opposite 
is time. By use of the dispersion relation 5»lt, the coefficient of the 
integral in equation 5.U may be written 
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which a simple examination demonstrates is always positive definite. 

The term under the integral therefore controls the sign of the contribu- 
tion at each pj^. Substitution of f^^ into this expression yields the 
relation 
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From symmetry condition B.U9* k/w is always positive definite for 
forward propagating waves. The pitch angle anisotropy that we are 
considering has Apj^ > Apii. It is concluded then that for + resonances 
expression 5*15 will always he less than zero. For minus resonances 
the sign of expression 5.15 will he controlled hy Jcv„/o). This tells us 
that plus resonances absorb energy from the waves when there is this sort 
of pitch angle distribution. The result is immediately suggestive of a 
mechanism for enhancing the high energy portion of the distribution. 

This is so because plus resonances only occur when (from equation 5*10) 
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When Y Aoes vary significantly from unity we are of course in 

j 

the relativistic regime. I 
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V.E, Scaling the Various Equations 

We have found that the numerous equations and relations that we 
have had to deal with are more easily analyzed in a scaled form. As 
we will he considering a numerical analysis it seemed appropriate to 
introduce the scaling at this point. We define the following unitless 
parameters 
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which form a consistent set that makes all our equations unitless. 
The dispersion relation becomes ; 
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The scaled linear growth rate is 
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and the resonance ccndition becomes: 
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Finally we write the scaled diffusion equation as: 
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in which the sum i rxins over the resonant values of K which satisfy 

the resonance condition for the particular values of p_, and p„. 

D, 


The 
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operator is defined hy 
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Dj^p X 3p„ X ^3pj^ 
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V.P. The Effects of Low Density 


As has been demonstrated by many authors, cyclotron \raves can be 
quite unstable in th'^ presence of a flat distribution. A significant 
amount of energy may be transferred from the particle distribution to 
the waves. We will demonstrate that at low post-substorm densities out- 
side the plasmapause that these waves will grow at the expense of the 
pitch angle anisotropy of the subrelativistic particles enabling the 
near relativistic particles undergoing plus resonances to feed on the 
increased wave energy. This will lead to an enhancement of the high 
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energy portion of the spectrum. 

Figure 5.2 is a diagram of (equation 5*22) for fixed values of 

Q 

as a function of x at a value of = ,1/cm (figure 5*2a) and 

Q 

= ,2/ cm (figure 5*213). The fixed values of p, increase upward. 

Plus resonances occur above the zero line of p« and minus resonance 
occur below '^he zero line. For a given frequency x the resonant 
parallel momentum, p,t, for a g’ven perpendicular momentum is a higher 
value for positive resonance for the n^ = ,1/cm case. The entire 
form of the resonant curves show that the positive values of p„ achieve 
greater values for the lower density case. The higher p,i for fixed 
then the higher is the energy. The conclusion is then that lower 
density tends to enhance the relativistic nature of the interaction. 

It pushes it to higher energy. 

V.G. Energy Diffusion and the Belativistic Case 

It is generally argued that diffusion in pitch angle is far more 
important for cyclotron resonant interactions (Kennel, I 966 ; Kennel, 

1969 ) than is energy diffusion. We would like to doaonstrate that energy 
diffusion can be quite significant in the regime we are considering. 

To facilitate this demonstration we will transform equation 5.23 to 
energy pitch angle coordinates. We take 
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The operator ^ — in energy pitch angle variables is expressed 
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which we write in mixed variables as: 


:— = 3j, a . /cos^a K X a 
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The term — is the inverse of the scaled wave phase velocity. For x 

less then ,1 the phase velocity is very small. Through most of the 

region of interest here it is of the order of 1/2 and fairly slowly 

varying. has a maximum value of 1, but as the energy, Y» 

increases the term decreases. Since most of the values are of 

yx 

order 1 , the first term is as important as the second. The first is 
energy diffusion, the second is pitch angle diffusion. The importance 
of the first term is emphasized however, at high energy and steep 
pitch angle, ^^hen this is coupled with the fact that plus resonances 
absorb energy in wave-particle interaction we see that we may expect 
to see important energy diffusion taking place for relativistic particles 
at steep pitch angles. 


V.H. Numericall?/^ Modeling the Quasi-linear Diffusion Equations 




The equations B.50, B.l »5 and 5*6 form a closed set which describe 
to evolution in time of the distribution function, f(Y»ct,t) (or 
equivalently f (Pj^sPti,t)) and the spectral energy density Ev('t)j of 
cyclotron waves. It is probably impossible to determine solutions to 
these coupled integrodiffcrcntial equations in a closed analytic form. 




rill ■iiTil^ri^ri.- 


iiATiBwai^ 
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We have however adopted a fairly standard procedure to determine the 
time "behavior of f(Y,a,t) numerically. f(Y»a,t) is established on a 
numerical mesh. Equation B.U5 is solved numerically for the growth 
rate at wave number K. An initial spectral energy density is assumed 
and equation 5*6 is solved numerically to deteimine after a 

small time increment At. With the new value of Ej^(t) equation B.50 
is evaluated to determine df/dt at each point on the mesh. The new 
f(YsCtjt + At) is determined by 

f(Y,a,t + At) ** f(Y^a,t) + ~ At 5.27 

With the new f(YsCe,t) the process is repeated. This procedure 
can he continued as long as desired. 

As formulated above the problom requires two dimensions and we 
have been unable to reduce it in any way to an equivalent one dimensional 
problem. The operator in equation (B.50) has numerous second derivatives 
inherent in it. Because of these second derivatives and the narrow 
range of instability a fine mesh was required. We were forced to smooth 
the first derivative in calculating second derivatives to remove obvious 
non-physical jagged behavior from the calculated second derivatives. 

The distribution function was defined over a 91 100 (9100) point mesh 

and the wave spectra was defined over 200 points. The dispersion relaticn 
was also defined over a 200 point set and logarithmic interpolation used 
to determine values between the mesh point. 

The initial distribution functions are chai*acterized by the parame- 
ters n^, Apj^ and Ap„. These values may bo varied to establish 
different initial growth rates and the location of the instability in 
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frequency. The waves must grow up against the convective loss to a 

level at which they affect the distribution function sufficiently before 

any change occurs in the initial distribution function. Because of 

the very complex nature of the numerical calculations a very large amount 

« 

of computation time is required. We have therefore only determined the 
value of ^ up to that time when ~ at any value of y and a 
becomes 1% of the initial value of f(Y>ct). Typically we have taken 
the ratio of backgrovmd particle density to effective superthermal par- 
ticle density to be 10:1. In figures 5.3a and 5«3b we show the ratio of 
change of the distribution function at the 1% point for the initial 
conditions 

o 

n = ,1/cm 
o 


( 1 ) 

Q 

ce 


.6 


Apit 


1.25 


5.28 


and Ap, chosen to give an equivalent perpendicular ten^erature of 
130 kev. To better represent the conditions in the magnetosphere we 
have established and maintained a three degree loss cone. All particles 
that appear in the loss cone are convected out of the interaction region. 
For clarity we have separated the region of negative change in the 
distribution {figure 5.^a) from the region of positive change in the 
distribution (figure 5.^tb). The results very clearly denonstrate the 
validity of our conclusions based upon the linear analysis. The region 
of greatest negative change in the distribution runs in an arch from 
about 100 kov and 80 degrees pitch angle up to 300 kev and 75 degrees pitch 
angle. The waves are undergoing negative resonance with those particles 






PIGUEE - RATE OF CHANGE OF DISTRIBUTION WITH COLD DENSITY 0.1 AND 
ANISOTROPY 1.25. RATIO OF ELECTRON FREQUENCIES 0.6. 
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and gaining energy from them due to their pitch angle anisotropy. The 
region of greatest positive growth (Figure 5.^h) runs in an arch from 3 
300 kev and 80 degrees pitch single up to 700 kev and 90 degree pitch angle. 
The waves are undergoing positive resonances with these particles and 
losing energy causing this region of the particle distribution to grow. 

As there cxe initially the least number of particles at high energy, the 
severest relative changes in the distribution will occur at hi^ energy. 
The inverse growth rates minus the convective time loss yield an effec- 
tive growth rate time for the instability. For the initial conditions 
5.28 this time is on the order of 1 second. 

The numerical results are in very good qualitative agi'eement with 
the experimental observations described in Chapter III. The relativistic 
particles are enhanced by the lower energy particles which should 
precipitate into the loss cone in a time like one over the growth rate 
(here about 1 second). We have of course varied the initial conditi.ons 
to guarantee that the resiilts would closely approximate to observations. 

A value of n^^ of .01/aar (.1 of n^) yields the right magnitude for the 
observed integral fluxes. The fact that a set of initial conditions 
for this problem could be chosen that approximate the observations maybe 
quite fortuitous in that we have only attempted a very simplistic 
(though mathematically complex) model interaction. It is really the 
qualitative featiires which may carry over into a more complete model 
incorporating the complexities of drifting trapped particles with 
possible sources and sinks. The spectral shape of the precipitated 
particles as well as the time behavior will depend strongly on the rest 
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of the model. The point ve have made here is that cyclotron resonant 
interactions vdll lead to the enhancanent of the relativistic electron 
population at times of low electron density. 

V.I. Conclusion 

We have not examined the precipitation of the relativistic electrons 
directly, hut only inferred that they should precipitate. A complete 
analysis should include a consideration of the electrostatic loss cone 
instability (F.VI. Perkins, 1959) as well as possible parasitic inter- 
action with background ion-cyclotron turbulence. It is possible to 
conceive of the foregoing theory as complementary to that presented by 
R. M. Thorne (1974). Thorne proposed a mechanism for producing REPs 
in which drifting protons injected on the nightside of the earth during 
a magnetic substorm would enter a socalled detached plasma region and 
would undergo strong pitch angle diffusion amplifying ion-cyclotron 
waves. The ' on-cyclotron waves would in turn interact with very ener- 
getic electrons causing them to precipitate (move into the loss cone). 

The mechanism of Thorne (1974) does not explain the production of the 
energetic elect:'"»ns but rather assumes they are present in the trapped 
particle flux. The theory does imply strong proton precipitation in 
association with the electron precipitation. Our theory explains the 
production of energetic electrons and does not require that there be 
accompanying protons (this agrees with our observations). In addition, 
our theory strongly couples the relativistic electrons to the precipitation 
of the lower energy electrons which also agrees very well with our 
observations. The gro\rth of the relativistic electron population through 
the anomalous doppler interaction on the other hand docs provide an 


input of relativistic electrons for the theory of Thorne ( 197 ^) and in 
that* sense the two theories complaient each other. 

There are two features of the REP theory which stand out as 
obvious indicators of its possible validity. The first is the association 
of very low cold plasma density at the equator with an ongoing REP 

< Sijjg) and the second a fairly narrow hand of strong whistler 
turbulence again near the equator. In ooir numerical model the magnetic 
fields in the turbulence reached a level of about 10 my over a band 
about 500 Ha wide before there was any appreciable change in the high 
energy portion of the distribution. Both of these measurements would 
have to be made on board a satellite. As there are presently and will 
be in the future satellites instrumented to make such measurements we 
may look forward to such measurements being made during an REP. 








APPMDIX A-1 - MODELING THE X-RAY PRODUCTION IN THE SHIELDING 

The experimental data on x-ray production to be foxmd in the litera- 
ture is displayed in many different forms and not all of it was directly 
useable for otir purposes. We have rescaled data taJcen from Berger and 
Seltzer (I968) and Rester, Dance and Davidson (1970), Figure A.l is a 
composite of this data for copper or materials with z values very close 
to copper- The x-ray spectra emergent between 0 and 10 degrees for 
electrons of various energies normally incident on targets of thickness 
R (the particle range as defined in section 2.D) is plotted in units of 
photons/(ker x ster x electrons) for electron energies between 200 Kev 
and 2.0 Mev. For purposes of this calculption we extrapolated those 
energies not actually covered by the data. 

Given the range R for various energy electrons the graphs in 
Figure A.l yield the x-ray spectra at the depth R , at which the 
electron cones to rest. The x-rays are then propagated the remaining 

distance through the brass, d - R , where d is the thickness of the 

o 

brass. The spectra is attemoated according to the law 

-Vp(d-R ) 

I(E,E) - I^(E,E^)e ® ° (A.l) 

with I^(E,E^) the intensity of x-rays of energy E at R^ and 
E is the energy of the original electrons producing the x-rays. The 
coefficients Uj, are tabulated in Liepunskii et al (1965) for numerous 
materials. The x-rays are attenuated again by the same law through the 
lead. For each incident electron energy there is an x-ray spectriua 
S(E,E) emergent from the lead. The x-ray spectrum is then folded with 
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The spec trim seen by the silicon detector will then be determined by the 
relation 


S (E) = S(E)(l-e 



(A. 3 ) 


where li„ is the appropriate coefficient for silicon also tabulated in 

Hi 

Leipxmskii et al (1965) and 6 is the detector thickness. 

The results of our mpdel calcoilations are described in section II . 2 
of the text. 







APPMDIX A -2 - EtlERGY SCALING PARAMETER, E 

r 

Many authors have scaled the transmission properties of thin foils 

for electrons in terms of an energy parameter, normally called E or 

r 

E , which makes the transmission curve nearly independent of the foil 
c 

thickness. This i^e of scaling is used through much of the literature 

on ionizing radiation and its various properties. Kanter (1961) credits 

Weaker (I 94 l) with being the first experimenter to notice this scaling 

property. Using Wecker's definition, Ranter's E is determined by 

c 

'drawing the tangent to the transmission curve at the point of steepest 
slope and taking the tangent’s intercept with the energy axis as the 
value of E^ . is supposed to be that energy at which transmission 

just begins. 

Berger and Seltzer in their various articles (see Bibliography) 
report the results of con5)uter simulation for which the above definition 
would be a bit awkward. They have adopted the scaling energy E^ 

(range energy) which is defined as that incident energy which yields a 
total electron path length, r , in using the continuous slowing down 
approximation (CSA). The values of E and E are very close with 

E generally being a bit smaller than E . We have chosen to use E 

C 37 * 

for all of our pxrrposes since it is an easily obtained CLuantity, The 

data we have taken from Kanter {1961) (Figure l.l6a) has been rescaled 

from E to E . 
c r 




APPBIDIX A-3 - AKGULAI? RESPOIISE OF THE COLLIMATOR 


For monodirectional flux of pEirticles of uniform intensity over the 
aperttire of the collimator it is simple to derive the response from geo- 
metrical consideration than from a vector analysis. We consider the flux 
as making a direction a with the normal to the plane of the back aper- 
ture (and front aperture) of the collimator, (see Figure A-2). The 
front aperture will either obstriict the back aperture, partially shade it 
or have it completely exposed. The exposed area of the back of the 
co]J.imator times the cosa will constitute the response of the colli- 
mator. For angles greater than 11.35 degrees the back aperture is 
con^jletely shaded and the response is zero. For angles less than 3.60 
degrees the back aperture is completely exposed and the response is 
sinply cosa (effectively 1). For angles between 3.6 degrees and 11.35 
degrees by simply projecting the shadow of the front aperture on the 
back apertxire it may be determined that 



- h sin0 < X. 

X 


6min < 9 < 6 

P 







I 


n 


: ) 


1 




.. ^ 

- 1 

l> WriM I MiiMM<rtif>ri>Mf;-i’ 1^11‘iift p < 1 . i 
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APPErmiX B - DERIVATION OF A FULLY RELATIVISTIC QUASILINEAR 
DIFFUSION EQUATION FOR CYCLOTRON RESONANT INTERACTIONS 


We followed the quasilinear fonnalxsm outlined in Krall and 
Trivelpiece (Krall and Trivelpiece, 1973) • We will assume that there is 
an interaction region centered at the equator of sufficient uniformity 


and size to define a uniform static magnetic field and a specially 


smooth distribution function 


~ V j ^ xf^(xj^,t) 


(B.l) 


so that perturbations satisfy the relation 


f^(x,p,t) = 


(B.2) 


and therefore 


1 1 d\f^^(x,p,t) = o 


(B,2a) 


We will make the same assumption for all perturbed field quantities 




d^x B 


(B.3) 


= 11 


E 


with 




i 


We substitute these values in. to the Vlasov equation and esqiand so that 



We take the spacial average and use condition B.2.a 






(B.5) 


Equation B.5 tells us that the time derivative of the spacial uniform 
distribution function is proportional to a product of perturbed quantities 
and is itself therefore a second order correction. The quasi~linear 
approximation consists in solving the first order equations for per- 
turbed quantities and substituting these into equation B.5 to determine 
df 

— rr~ . Dropping second order terms then in equation B.U we have the 

CLv 

standard perturbed Vlasov equation 

(B.6) 



As we wish to retain all relativistic terms we solve this completely 
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here using the method of characteristics which cw he found in most 
standard textbooks on Plasma Physics (see Krall and Trivelpiece). 
One recognizes that equation B.6 may be written 


df. 


al 


dt 






(B.7) 


where the operation “ is evaluated along the trajectories of unper- 

CLti 

turbed particles in the tiniform magnetic field. We integrate both sides 
of B.7 from - to t and assume that the perturbed qiiantities were 
zero at - “ . We then have that 


f / v»xB!,\ 3f 

where we must remenflier that E, and 3, are functions of xlv' and p the 

XI — — 

parameters of the unperturbed trajectories. We now choose and ^ 

in such a manner that 


x‘ (t' = t) = X 

V'(tf = t)=v (B.9) 

p‘(t'=t) = p 

VS ^ 

We take k along the magnetic field and (j) as the azimuthal direction 

A ^ 

V = V„ Xl ^ j 

P = + Pi.? I 

• ^ 
I 


about the field with 
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then 

z* = z + v« (t'~t) 

x'= X - sin(4>-u (t*-t)) + sinifi 
“ca ca 

y'= cos(4>-w (t'-t)) - cos4> (B.IO) 

ca “ca 

= PACOs(^-Ul^^(t'-t)) 

where (i) is the cyclotron frequency for the paa^icle of mass m such 
cot ct 

that 


q B 

W =1-2. 

ca vm c 
a 


(B.ll) 


— 1/2 

and Y = (l - “ 2 )” ' is the relativistic mass correction. 

c 

If we now replace all functions of x and x' in B.8 hy their 
Fourier integral representations such that 


f^(t.2£> 


-JL_ f 

(2«)^ > 




(B.12) 


and change the order of integration of k and ib eqtiation B.8 "becomes 


_ / 4 .\ ik*x 



8f 



ik*x* 

e 


dt' 


(E.13) 
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We repeat the procedtire over all fmctions of t and t* with 

* JC”“i 

and d vide the resulting equation hy e — — . Equation B.13 now 

■becomes 

(B.15) 
t 

We*“> = - 1 1 






Sf 


ao ik. (x’-x) -io3(t'-t) .,i 

Tr-“ e ~ ” e at 

3p 


It is convenient to let t* - t = j and x| - x = Ax in equation 
B,15» The orientation of the x and y axies is arbitrary so that we 
may let = k and k^. = o without any 7oss of generality. 

The term in the exponential on the ri^t side of B.15 is expanded 

(B.16) 

k * ^(t) - wt = (sii.>sin(<{)-a)^^T) } + (k„T^, - u)t 

‘“ca ^ 

To evaluate the integral the vector term must be expanded. Using 
Maxvrell's equations the perturbed quantity related to the 

perturbed quantity ^ (w) 

(w) 

B (w) (B.17) 

—Ik t«i 

Equation (B.l?) is substituted into (B.15) and the vector product 


expanded, yielding the result 
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(B.20 conii) 


ne*“® 


n=-<° 


Wiere is the Bessel fvir-ction of the first hind of order n , The 


three types of integrals are then written 


(B.21) 


0 

1 


dT e 


ibsin4> -ibsin(ij>-a) x) i(hjT 5 -tij) 


ca e 


C O o 


_ ^ihsin<j) 


£ f at J (b) e-^'* 

n=-^ 


f at e*’’^*”* ^-lisln(<ii-«^^T) ^i(lt.T. -u)T =. 


_co 


ie' 


ibsin^) 


00 

o 

Z_» f — (j 
.IrL J 2 


(b)-J ^,(b)) e-l"* ax 

n +1 


I ax cos($-«„x) 


ca 


UX> 




jl=^<n _ioa 


Where we have tolten b - 


Ui 

ca 
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Equation B.l8 may then he reduced to 






; -iE +E - -xE -E .. 

2 J 2 1 -0)^^ -w _ 


J 


with 


k V 3f 
_ g go 

'1 w 3p„ 


o, = -rr— + (1 - — ) 2^ 


3f 


3P, 


oo 

2 


Using the standard transformation for left and right handed circularly 
polari'ied waves we let 


and 


®L = 


E H* 

E - iE 



f (PiW) = o < 
olk^^’ 2 11 


-i<}> 




i(J> 




(B.25) 


(B.26) 


f^j^(^,u)) from equation B.26 now defines the first order perturbed 
quantity which is now placed back in equation (B.5) to determine the 
second order time derivative of the initial spacially averaged distri- 
bution function. Substituting the Fourier integrals on the right hand 
side of equation (B.5) we have 


i 

f 



V(2w)^ 



,3, ,3,1 
d kd k 


I 


, , , vxBik(w) 

(E^(w) + = ) 


{B.27) 
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^^caK* ik»x i}i‘X JS 
— 


(B.27 cont.) 


and ve identify 


1 f ik.x ih?*x ,3 

— S' e e dx = 

} 


dCk-f-k*) 


with the resxilt 


I AA /e (ttj) + ^ 

} V l%k^ ^ c 


^W“> 


(B.28) 


Since we require that he resil, then 


^alk 

Again we expand the vector produce replacing B^^^((o) in terms of 

E (w) . Using the relation 
"Ik 


E^ cos4> Ey si»4> ~ 2 ^ 


the vector product becomes 


^Ejj e-^^CA-lB) * Ej_ e^'^(A->'iB)j- 




with the operators 


0, 3p„ 


B = (l - ~ ~ 


{B.31 cont.) 


Evaluating the terms in B.31 with equation B.26 yields 


I Ej, 




k V "^W“ -to 

^ » ca 


(B.32) 


- h I (1 - ^)] 

I Ejj e-i*(A-iB) = =|a|[A * i (1 - Sj^)] 


k v -u) 
« » ca 




k,,v, “W -w 
u II ca . 


The parameter ^ was the arhitrary initial phase angle between Pj_ 
and the x-axis. This will be assumed to be a random quantity. We will 
integrate both sides of equation (B.28) then over <fi from o to Sir 


and divide both sides by 2 tt 


(B.33) 


df 
go _ 

dt “ Sir 


’ yj® (w)\ 3f* (w) 

%k<“) * 


and note that 


+2i«i> 


We now have 



We will use the standard definition of the spectrsil energy density 
(Krall and Trivelpiece, 1973) with 
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The dispersion relation which relates sa and k in equation B.37 must 
be determined before B.37 may be evalmted 

The dispersion relation is determined by standard techniques of 
linear plasma theoiy. Using Maxtrell's equations with the assumption of 
transverse waves it is easy to determine that 




Utt p 

ill) i 


I / ' 


(B.38) 


where is the density of the a-species. 


With V =* Vi. (e^^ Sp + e_ ) + ^ 


E = Cjj + Ej^ where e^^ and e^^ are the basis vectors for the 
right and left hand circular representations. Equating vectors component 
wise (B.38) becomes 




^3 i4» ^ 


(B.39) 




We will integrate over (}) is equations (B.39) in which 


dp = p^ dp^ d^) 


as in equation (B.32) cross terms of the form and drop 
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out and ve are left with the Landau dielectric for R and L mode 
cyclotron waves 




I Pi.dPiap„ 




w 3p^ 


L.^YJmX 1 _ 

V w y 3p^j -w 


The plasma freq.uency of species-a is given "by 




pa 


The standard approximation requires that so that D(kjW) may 

he expanded in a power series in giving the value of Uj as 


“l 9D. 


We have taken D = D_ + iD_ . The imaginary part of D(k»(o) comes from 
the p integration in equation (B.ltO) where the relation 


0) '♦o 
1 


im f dp r — 


f kv, +to -u) 
■» “ ca 


iTiym^^F(^) 


P«=Pr 


f w- \ 

with p as the Cauchi principal value and = yma determines 

D_. . 
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Using relations B.UO, B.Ul, and B.It2 ve detenaine 

(B.U3) 



At this point we consider that we have three plasma components, a 
relatively cold electron and proton gas (at ahont 2000°k) and a super 
thermal electron component. We will assume thau there is over all 
charge neutrality with 


and 


n n 
ec pc 


n 

0 


"eH = “l “c 


(B.43) 


where n , n and n „ are the respective densities of the cold 
ec pc eH 

electron, cold proton and super thermal electron componc-its. Conditions 

B.U3 allow us to take the real part of the dielectric D„(k,w) strictly 

r\ 

from the cold plasma components and the imaginary part of the dielectric 
Dj strictly from the super thermal component, (Montgomery and Tidman, 
196U), We arrive finally then at the relations 


0 )^. 1,2 2 

ce Cl 


(B.liU) 


eB 

and with 

cd m c 




and 




^-kL * ®kR 
^-kR ” Sl 

(Kt€i 11 and Trivelpiece, 1973), we see that the terms mder the integral 
cancel except for their imaginar?^ parts. 

I 

To this point we have retained the LH-mode denominator "because as 
we see it provides the term necessary to cancel the real term under the 
integral in equation (B.li8). This guarantees that the term on the ri^t 
hand side is completely real, ^•Jhen we determine the imaginary part of 
with the help of relation (B.U2) we obtain four terms; two terms 
for RH-mode and two terms for LH~mode. At this point we nay drop the 
IiH-mode since it makes no further contribution. Equation (B.U8) now may 
be written (again k = k) 

0 


(3.50 
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TJie first term (which is lander the integral) is refered to as the 
non-resonant diffusion tern and the second term (which is not under the 
integral) is refered to as the resonant diffusion term. 

We use eq,uations (B. 50 ) and (B,45) as the starting point for our 
discussion in the text. 


I 


! 


I 


— } 



i 


j 

I 

! i- SEIECTED BIBLIOGRAPHY 


Akasofu, S., Polar and Maf?raetospheric Substorms, Springer-Verlag, Hew 
York Inc. , New York,' 1968, 

Ansari, Z. A., A peciiliar type of daytime absolution in the auroral 
zone, J. Geouhys. Res. , 70, SUTs 1965. 

Bailey, D. K. and M. A. Pomerantz, Relativistic electron precipitation 
into the mesophere at subauroral latitudes, J. Geophys. Res ., TO, 

5823, 1965. 

Bailey, D. K. , Some quantitative aspects of electron precipitation in 
and near the auroral zone. Reviews of Geouh ., 6, 289, 1968. 

Bendat, J. S. and A. G. Piersol, Heasxnrement and Analysis of Random Data, 
John-Wiley and Sons, Inc., New York, 1966. 

Berger, M. J. , and S. M. Seltzwe, Studies in penetration of charged 
particles in matter. Nuclear Sci. Series, Report #39, Commission on 
Nuclear Science, National Academy of Science, HRC Commission on 
Nuclear Science pub, 1133, 1964. 

Berger, M. J. , and S. M. Seltzer, Penetration of electrons and associated 
Bremsstrablung through aluminum targets. Protection Against Suace 
RadlatioE, NASA. SP~l69, 2o5, 1968. 

Berger, M. J. and S. M. Seltzer, Bremsstrahlung in the atmosphere, 
Atmospheric and Terrestial Physics , 34, 85, 1972. 

Berger, M. J. , S. M. Seltzer, and K, Maeda, Energy deposition by auroral 
electrons in the atmosphere, J. Atmospheric and Terrestial Physics, 

^ 1015, 1970. 

Berger, M, J. , S. M. Seltzer, and K. Maeda, Some new results cn electron 
transport in the atmosphere, to be published in J. A t mospheric and 
Terrestial Physics, 1975. 

Berger, M, J., S. M. Seltzer, S. E. Chappell, J, C. Humphreys, and J. H. 
Motz, Table of Response Functions for Silicon Electron Detectors , 

UBS Technical Note 489, 1969. 


Bethe, H. , Theorie des durchgangs schneller korpuskular-strachlen durch 
materie, Annelan der Physik, _5, 325, 1930. 


6 - 



Bogott, F. H. , and F, S. Mozer, Drifting energetic particle bunches 
observed on ATS 5, J. Geonhys . Res . , 75 , 1825, 1974. 

Brice, H., Harnessing the energy in the radiation belt§, J. Geophys . 

Res ., ;j6 , 4698, 1971. 

Brice, N., and C. Lucas, Influence of magnetospheric convection and polar 
wind on loss of electrons from the outer radiation belts. 


j 

■■'i 






182 


J. Geophys, Res ., 76, 900, 1971. 

Buechner, W. W. , E, J, Vaa de Graaff, E, A. Burrill, and A. Sperduto, 
Thick-target x-ray production from 1250 to 2350 kilovolts. Physical 
Eeview, jU, 1348, 1948. 

Carpenter, D. L. , Whistler studies of plasmapause ir. the magnetosphere, 

J, Geophys, Res. , 71, 693, I966. 

Carpenter, D. L,, and R. L. Smith, vnaistler measuretients of electron 
density in the magnetosphere. Reviews Geophys . , 2^, 4l5, 1964. 

Chappel, C. R., K. K. Harris, and G. ¥. Sha3rp, A study of the influence 
of magnetic activity on the location of the plasmapause as meastired 
hy OGO 5» J. Geophys . Res . , 75, 50, 1970. 

Chappel, C. R., K. K. Harris, and G_. W. Sharp, The dayside of the 
plasmapause, J, Geouhys. Res. , 7*6, 7632, 1971* 

Cocke, ¥. J., and J. H. Cornwall, Theoretical simulation of micro- 
pulsations, J. Geophys. Res . , 2843, 1967- 

Cornwall, J. ¥. , P. V. Coroniti, and R. M. Thome, Turhulent loss of 
ring current protons, J. Geophys. Res ., 75, 4699, 1970. 

Coroniti, F. V., and C. P. Kennel, Electron precipitation ptilsations, 

J. Geophys. Res ., 75 , 1279, 1970. 

Coroniti, E. V., R. ¥. Fredricks and R. "Sf/hite, Instability of ring 
current protons beyond the plasmapause during injection events, J^. 
Geophys . Res . , 77, 6243, 1972. 

Davidson, R. C., and D. A. Hammer, Nonequilibrium energy constants 

associated with large-aa^litude electron whistlers. Physics of Flxiids, 

15, 1282, 1972. 

Dearnaley, G. , and D. L. Northrop, Semiconductor Counters for Nuclear 
Radiations, John Wiley & Sons, Inc., New York, I966. 

Dessler, A. J., Solar wind interactions and the magnetosphere. Physics 
of the Magnetosphere , Ed. R, L, Carovillano, J. F. McClay, and H. R. 
Eadoski, D. Reidel Publishing Company, Dordrect, Holland, 1968. 

Drummond, N. E. , Quasilinear theory of plasma turbulence. Plasma Physics, 
International Atomic Energy Agency, Vienna, I965. 

Drummond, N. W. , and D. Pines, Non-linear stability of plasma oscillations. 
Nuclear^ Fusion, Surnl, Part 3 , 1049, 1962. 

Engle, R. D. Non-linear stability of the extraordinary wave in the 
plasma, Phy.^,ics of Fluids , 8, 939* 1965 - - 

Felds toin, Y, I., Polar auroras, polar substonas, and their relation- 
ships with the djTiamics of the magneto sphere, Proceedings of the 


International Syinpositon on the Physics of the Magnetosphere ^ Ed. D. J, 
Williams and G. D. Mead, American Geophysical Unions 1968. 

Gendrin, R. J. , Pitcli angle diffusion of low energy protons due to 
gyroresonant wave interactions, J. Atmospheric and Terrestial 
Physics, 30, 1313, 1968. 

Goldstein, H., Classical Mechanics, Addison-Wesley, Wew York, 1950. 

Hardy, A. C., M. D. Lopez, and T. T, White, A parametric technique of 
computing primary electron dose, Trans. American ITuclear Soc. , 10, 

383, 1967. 

Hass, G., W. R. Hunter, and R. Tousey, Reflectance of evaporated 
Aluminum in the vacuum tiltraviolet , J, Optical Seo. American, 46, 

1009, 1056. 

Hass, G,, and J. E, Waylonis, Optical constants and reflectance and 
transmittance of evaporated aluminum in the visible and ultraviolet, 

J. Optical Soc. America, 51 , 719, 19^1 • 

Hillas, A. M. , Cosmic Rays, Pergaraon Press, Oxford, 1972* 

Kanter, H. , Electron scattering by thin foils for energies below 10 kev. 
Physical Review, 121, 46l, I96I. 

Kennel, C. F. , Low frequency whistler mode. Physics of Fluids, £., 219O 

1966. 

Kennel, C. F. , Consequences of a magnet ospheric plasma. Reviews of 
Geophys . , 7, 379, 1969. 

Kennel, C. F., and R. D, Engleraaiin, Velocity space diffusion from weak 
plasma turbulence in a magnetic field, I^ysics of Fluids, 9, 2377, 
1966. 

Kennel, 0. F. , and H. E. Petsch^, Limit on stably trapped particle 
fluxes, J. Geophys. Res. , 71, 1, 1966, 

Kobetich, E. J. , and R. Ifetz, Energy deposition by electron beams and 
delta-rays. Physical Review, I70, 170, I968. 

Kock, H. M. , and J. W. Motz, Bremsstrahlmg cross-section formulas and 
related data. Reviews of Modem Physics, 31, 920, 1947. 

Krall, N. W., and A. W. Trfvelpiece, Principles of Plasma Physics 
McGraw Hill Book Co., New York, 1973. 

Leipunskii, 0. I. B, V. Novozhilov, and V. H. Sakharov, The Propagatio n 
of Gainma Quanta in Matter , Vol. 6, Pergamon Press, 1965. 

Lerche, I., Space charge waves and isotropy of the Galactic cosmic ray 
plasma. Physics of Fluids , 10, I07I, 1967* 








184 


Lezniak, T. W, , and J. R. Winckler, Structure of the magnetopause at 
6,6 R„ in terms of 50 to 150 kev electrons, J. Geophys, Res., 73 
5733, 1968. 

Lezniak, T, W. , and J, E. Winckler, Experimental study of magnet ospheric 
motions and the acceleration of energetic electron during suhstorms, 

J, Geophys. Res ., 75, 7075, 1970. 

Liemohn, H. , VLF and ULF whistler anplification, J, Geophys. Res . , 72, 

39, 1967 . 

Little, C, G. and H. Lembaeh, Some measurement of high-latitude iono- 
spheric absorption using extraterrestial radio waves, Proc. I. R. E . , 

46» 334 , 1958 . 

Lyons, L., and R. M, Thorne, The magnet ospheric reflection of whistlers. 
Planetary and Space Science , I 8 , 1753, 1970. 

Lyons. L. , and R. M. Thorne, Parasitic pitch angle diffusion of radiation 
belt particles by ion cyclotron waves, J. Geonhys, Res . , 77, 5608 , 
1972. 

Lyons, L., R. M. Thome and C. F, Kennel, Electron pitch-angle 
diffusion driven by oblique whistler-mode tTirbulence, J. Plasma 
Physics , 589 , 1971 . 

Lyons, L., R, M. Thome and C. F. Kennel, Pitch-angle diffusion of 
radiation belt electrons within the plasmasphere, J. Geophys. Res. , 

n, 3455 , 1972 . 

Meada, K. , Diffusion of low energy electrons in the atmosphere, J,. 
Atmospheric and Terras tial Physics, 27, 259, 1965. 

Matthews, D, L., and D. J. Sisons, Observations of relativistic electron 
precipitation at L=6, J. Geophys. Res., 78 , 7539, 1973. 

Miller, W. , J. W. Motz, and C. Cialella, Thick-target bremsstrahlung 
spectra for 1.00, 1.25, and 1. 40 Mev electrons. Physical Review, 96 . 
1344, 1954 . 

Montogomeiy, D. , and D. A. Tidman, Plasma Kinetic Theory, McGraw-Hill, 
Hew York, 1964. 

Mozer, F. S., Particle-flxix limits in synchronous orbit, J, Geophys. 

Res. , XL> 2401, 1972 . 

Ness, N. F., The geomagnetic tail, Marnato spheric Fhysics; Froc. of the 
International S^^nposi^lm on the Pliyslcs of the Magnetosphere, Ed. 

D. J. Williams and G. D. Mead, American Geophysical Union, I 968 . 

Northcliffe, L. C., and R. F. Schilling, Range and stopping power tables 
for heavy ions, Nuclear Pain Tables A7, 233, 1970. 


185 


Ortner, J, , and W. Riedler, A smooth tyjie of cosmic noise absorption, 
Nature, 20l*, H8l, 1964. 

Ossakow, S. L. , E. Ott and I. Haber, Theory and computer simulation of 
whistler turbulence and velocity space diffusion in the magnetospheric 
plasma, J. Geophys. Res. , T8, 2945, 1973. 

Parks, G. K. , and J. R. Winckler, Acceleration of energetic electrons 
observed at the synchronous altitude during magnetospheric sub- 
storms, J. Geophys, Res ., 73, 5786, I968. 

Pfitzer, K. A,, and J. R, Vinckler, Experimental observation of large 
addition to the electron inner radiation belt after a solar flare 
event, J. Geobhys. Res. , 73, 5792, 1968. 

Rester, D. H., W. E. Dance and J. H. Derrickson, Thick target 

bremsstrahlung produced in electron bombardra^t of targets of Be, Sn, 
and Au in the energy range 0.2-2. 8 Mev, Journal of Applied Physics, 

4i, 568, 1970. 

Roberts, C. S., Pitch-angle diffusion of electrons in the magnetosphere. 
Reviews of Geophysics and Space Physics , 7^, 305, 1969- 

Roberts, C. S., and M. Schulz, Bounce resonant scattering of particles 
trapped in the Earth *s magnetospheric field, J, Geophys. Res . , 73 , 

7361, 1968. 

Rosenberg, T. J,, L, J. Lanzerotti, D. K, Bailey, and J. D. Pierson, 
Energy spectra in relativistic electron precipitation events, 

Joxinral of Atmospheric and Terrestial Physics , 34 , 1977, 1972. 

Rossi, B., and S. Olbert, Introduction to the Physics of Space , McGraw- 
Hill Book CoB^any, T970. 

Roux, A., and J. Solomon, Self-consistent solution of the quasi-linear 
theoiy. Application to the spectral shape of the VLF waves in the 
magnetosphere. Journal of Atmospheric and Terrestial Physics, 33, 

1457, 1971. 

Sen, H. K. , and A. A, Iflyller, Generalization of the Appleton-Hartue 
magneto-ionic formula, Phys. Rev. Lett ., j+, 355, i960. 

Stassinopoulos, E. G, , World Maps and Constant B, L, and FLUX Contours , 
NASA-SP-3054, NASA, Greenbelt, Md., I968. 

Taylor, J. M. , Semiconductor Particle Detectors , Butterworth Inc., 
Washington, 1963. 

Thome, R, M., Unducted whistler evidence for secondary peak in the 
electron energy spectrvim near 10 kev, J. Geophys, Res,,- 73, 4895, 

1968. 

Thorne, R. M, , A possible cause of dayside relativistic electron 




186 


precipitation events. Journal of Atmospheric and Terrestlal 
Physics , 36 , 635 » 197^ 

Kiome, R. M., and C. F, Kennel, Relativistic electron precipitation 
during magnetic storm main phase, J. Geophys. Res ., 76 , 1971* 

53iorne, R. M, , and C. F, Kennel, Locating the magneto spheric ring 
current, Connnents Astronhys, and Space Science ,. 3^ 115, 1971. 

Thome, R, M. , E. J. Smith, R. K. Burton, and R. E. Holzer, Banded hiss 
em' 'sion from the plasmapause, J. Geonh^/s. Res., 78 , I 58 I, 1973. 

Tvershoi, B. A. , Theory of turbulent acceleration of charged particles 
in a plasma, JETP, 26, 1032, I 963 . 

Van Allen, J. A. , Carged particles in the Magnetosphere, Magnetospheric 
Physics t Proc. of the International, Symposium on the Physics of the 
Magnetosphere , Ed, D. J, Williams and G. D, Mead, American Geo- 
physical Union, I 968 . 

Wecker, F. , Neue messungen der absorption ruckdiffuslon und 

sekundarstzahlenegTxag mittelschneller kathodenstrahlen in aluminum 
und gold, Annaleon der Pfaysik, Uo, 405, 19^1 . 

Whittaker, E. T. , and G. H, Watson, A Course of Modem Analysis, 
Cambridge Ifaiversity Press, 1965* 

Wiley, R. E., and F. D. Barish, Plots of Geomagnetic Field Geometry, 
Technical Report Ho, AFl-IL-TR-69-144, Air Force Weapons Laboratory, 
Albuquerque, Hew Mexico, I 966 . 

Williams, D. J. and L. R, Lyons, The proton ring current and its 

interaction with the. plasmapavise: storm recovery phase, J, Geophys . 

Res ,, 79 , 4195, 1974. 

Williams D. J. , and L. R. Lyons, Further aspects of the proton ring 
current interaction with the plasmapause: main and recovery phases, 

J, Geophys, Res ., (to be published 1975). 

Williamson, C., and J. P, Bougot, Tables of range and rate of energy 
loss of charged particles of energy.. 5 to 150 Mev, Rappt, CEA 2189, 
Centre Etudes Hucl., Saclay, France, I 962 , 

Yakimenko, V. L., Absorption of waves in a plasma (quasi-linear 
approximation), JK'l'P , 17 » 1032, 19^3* 


.. * 


r mmrtii II< um I Ti 






VITA 


Name: David Jay Simons 

Permanent address: 85 Joya Loop, Los Alamos, New Mexico, 875^^ 


Degree and date to be conferred; Ph.D. , 1975 
Date of birth: 

Place of birth: 

Secondary education: Forest Hills KL^ School, New York, N. Y., 

U. S, Grant High School, Los Angeles, Calif., I 96 I. 


Collegiate institutions attended 

Itoiversity of California, Berkely 
Ifeiversity of Maryland, College Park 


Dates 


Degree Date of Degree 


1961-1966 AB 
1969-1975 Ph.D. 


1966 

1975 


Major: Astrophysics and Geophysics 

Minor: General Physics 

Publications and Presentations to Learned Societies; 


Mhtthews, D. L. and D. J. Simons, Observations of relativistic electron 
precipitation at L=6, J. Geophys. Res., 78 , 7539* 1976. 

Simons, D. J. and D. L. Matthews, Detailed Measurement of nergetic 

electrons in a REP event, EOS, 5^ Transactions, American Geophysical 
Union, U 2 I, 1973. 

Pongratz, M. B., K. Papadopoulos , and D. J. Simons, Observations of the 
bump-in tail instability in auroral electrons, EOS , Transactions, 
American Geopl^sical Union, 54, 421, 1973. 

Simons, D. J., and D. L. Matthews, Production of relativistic electrons 
on L=6 during times of post substorm electron density, EOS, 56 , 
Transactions, American Geophysical Union, 1174, 1974. 

Pongratz, M. B., J. Wolcott, C. Keady, and D. J. Simons, Aircraft ob- 
servations of equatorial mightglow emissions of atomic o^gen 6300 A 
and 5577 EOS, 56 , Transactions, American Geophysical Union, 1152, 

1974. 

Positions held: Jan. 1963 to April 1967» Research Data Analyst, Lain'once 

Radiation Laboratory, University of California, Berkeley, 
California. 

April 1967 to April 1969 * United State Peace Corps 
Volmteer, U. S. Peace Corps, Lagos, Nigeria. 

88ECEDING PAGE BLANK NOT PILteD 




Positions held; September 19o9 to Jan. 1971» Graduate Teaching 
Assistant, Department of Physics, University of 
Maryland, College Park, Maryland. 

January 1971 "to September 197^ s Graduate Research 
Assistant, Institute for Fluid Dynamics and Applied 
Mathematics, University of Maryland, College Park 
Maryland. 

September 197^ to present. Senior Scientist, EG & G 
Inc., Los Alamos, Hew Mexico. 



